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FOREWORD 


This  report  is  the  second  in  a  series;  the  first  report  is 
ESD-TDR-64-112,  "Engineering  Simulation  of  Powered 
Flight,"  by  R.  W.  Dix,  The  MITRE  Corporation,  July  1964. 


ABSTRACT 


This  report  is  the  second  of  a  series  describing  the  current  effort  towards 
establishing  a  workable  engineering  simulation  of  the  space-ground  environment. 
This  report  describes  the  simulation  of  ground  station  coverage  of  a  vehicle 
from  the  end  of  powered  flight  to  an  arbitrary  time  in  orbit.  Signal  strengths  of 
several  vehicle-ground  station  antenna  combinations,  as  well  as  geometrical 
coverage,  are  simulated.  A  spherical  earth  model  was  used  in  the  simulation, 
and  the  effects  of  atmospheric  drag  and  orbital  perturbations  were  neglected. 


REVIEW  AND  APPROVAL 


Publication  of  this  technical  documentary  report  does  not  constitute  Air  Force 
approval  of  the  reports  findings  or  conclusions.  It  is  published  only  for  the 
exchange  and  stimulation  of  ideas. 


/  /  JOHN  N.  VAN  DUSEN 
'  /  Major,  USAF 
!  /  ESLN 
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SECTION  I 


INTRODUCTION 


This  report  is  the  second  in  a  series  describing  the  current  effort  toward 
establishing  a  workable  engineering  simulation  of  the  space-ground  environment. 
The  activity  divides  into  three  phases: 

Phase  1;  Powered-flight  and  free-flight  trajectory  simulation. 

Phase  2:  Simulation  of  the  interactions  between  the  space  vehicle, 

other  vehicles,  and  the  ground  environment,  such  as  tracking  cover¬ 
age,  telemetry,  communications,  etc. 

Phase  2:  Engineering  utilization  of  the  above  simulations  in  the 

solution  of  a  particular  problem,  such  as  evaluating  or  optimizing 
the  communication  capacity  of  a  given  system,  or  evaluating  the 
ground  environment  in  relation  to  a  given  set  or  type  of  missions. 

The  powered-flight  trajectory  simulation  of  Phase  1  has  been  documented  in 
ESD— TDR-64-112.  The  free-flight  trajectory  simulation  of  Phase  1  and  most  of 
Phase  2  is  documented  in  this  report.  More  specifically,  the  program  described 
in  this  report  calculates  the  orbital  history  of  a  space  vehicle  from  the  point  of 
injection  to  a  selected  time  in  the  future,  the  tracking  coverage  of  that  vehicle 
by  an  arbitrary  number  of  ground  stations  from  launch  through  orbital  history, 
and  the  electromagnetic  signal  strengths  of  various  space-vehicle  ground-station 
antenna  configurations. 
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PROGRAMING  APPROACH 


The  engineering  simulation  of  the  space-ground  environment  was  developed 
as  a  series  of  programs,  each  of  which  performs  a  specific  engineering  function. 
This  approach  allows  considerable  flexibility  and  adaptability  to  new  situations. 

The  over -all  simulation  is  presently  separated  into  five  programs  which  are 
related  as  shown  in  Fig.  1.  One  program  (MAIN)  provides  executive-tj'pe  con¬ 
trol  and  acts  as  the  main  calling  routine.  The  other  four  programs  (TRAC, 

RUNG,  CALC,  DEBE)  perform  specific  engineering  simulations  and  are  called 
subroutines.  The  programs  are  described  in  Section  2  under  the  heading  General 
Description. 

TYPICAL  APPLICATION 

This  program  was  used  in  the  study  of  emergency  communication  vehicles. 
Various  trajectories  were  simulated,  and  the  trade-off  between  payload  weight 
(i.  e.  ,  transmitted  power)  and  the  number  and  distribution  of  the  ground  stations 
that  could  communicate  with  the  payload  were  evaluated.  It  was  possible  to 
quickly  arrive  at  a  trajectory  and  payload  weight  which  satisfied  the  operational 
requirements  of  the  systems.  The  program  was  also  used  to  evaluate  the 
available  tracking  coverage  for  the  flight-test  trajectory. 

ORGANIZATION  OF  THE  REPORT 

The  report  has  been  arranged  with  the  interests  of  the  reader  in  mind.  An 
effort  was  made  to  include  enough  information  in  Section  2  to  permit  someone 
with  a  general  knowledge  of  digital  simulation  to  use  the  program  without  further 
reading.  If  a  deeper  understanding  of  the  simulation  is  desired,  the  reader  will 
find  the  mathematical  derivations  corresponding  to  each  of  the  subprograms  in  the 
remaining  sections.  The  appendixes  contain  a  complete  set  of  Fortran  IV  listings. 
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I 


CASE  LOOP 


Fig.  1  Functional  Flow  Diagram 
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SECTION  II 


USER  INFORMATION 

The  program  described  in  this  report  assumes  a  spherical-earth  model; 
the  effects  of  earth  perturbations  and  atmospheric  drag  on  the  orbit  are  neglected. 
As  originally  conceived,  this  program  was  used  in  the  study  of  emergency  com¬ 
munication  rocket  systems  involving  high-loft,  suborbital  trajectories.  Since 
trajectory  times  were  relatively  short,  earth  perturbation  effects  could  be 
neglected;  considerable  programing  time  was  saved  by  assuming  a  spherical- 
earth  model.  For  orbital  trajectories  involving  longer  periods  of  time,  earth- 
perturbation  effects  cannot,  in  general,  be  neglected.  The  user  must  decide  if 
these  effects  are  important  in  his  particular  application. 

GENERAL  DESCRIPTION 

The  detailed  interrelation  of  the  five  programs  is  shown  in  the  flow  chart 
given  in  Fig.  2. 

Program  MAIN  acts  as  the  executive  or  main  calling  routine.  The  param¬ 
eters  defining  the  vehicles  position,  velocity,  and  attitude  at  burnout  of  the 
final  stage  of  the  booster  are  read  into  the  program.  Subroutine  TRAC  is  called. 
MAIN  also  provides  a  logical  Fortran  end  of  program  if  the  burnout  velocity  is 
negative.  This  provision  is  a  programing  formality  and  normally  is  never  used. 

A  normal  exit  occurs  when  the  input  data  have  all  been  read  in. 

Subroutine  TRAC  calculates  the  orbital  position  relative  to  axes  fixed  in  a 
spherical,  rotating  earth.  It  also  controls  the  calling  of  the  other  subroutines. 

If  station  coverage  is  required,  CALC  is  called;  if  station  data  is  not  desired,  the 
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Fig.  2  Functional  Flow  Diagram — Basic  Orbit  and  Sensor  Coverage  Program 
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time-varying  orbital  data  is  printed  out  in  TRAC.  The  basic  control  data  card 
for  this  option  and  other  options,  time-interval  controls,  and  integration  sequencing 
is  read  in  by  TRAC.  These  control  constants  are  defined  below  in  the  description 
of  input  data  cards. 

Subroutine  CALC  performs,  in  general  terms,  the  following  functions: 

(1)  Reads  tracking  station  geographic  and  antenna  data.  Up  to  150 
stations  are  permissible. 

(2)  Reads  vehicle  and  sensor  antenna  performance  data. 

(3)  Calculates  geometrical  relations  between  instantaneous  vehicle 
position  and  stations. 

(4)  Calculates  rates  of  change  of  the  above  relations  by  time  averages. 

(5)  Calculates  angles  defining  several  vehicle-sensor  antenna  axis 
orientations. 

(6)  Calls  DEBE  option. 

(7)  Prints  out  orbital  and  station  coverage  data. 

Subroutine  DEBE  calculates  antenna  performance  and  space  losses  for 
eight  combinations  of  vehicle  and  ground  antenna  configurations.  The  vehicle 
antenna  may  be  either  inertially  fixed  (spin  stabilized),  or  constantly  pointing 
towards  the  center  of  the  earth  (earth  centered),  or  it  may  be  tracking  the  ground 
stations  (vehicle  tracking).  The  ground  antenna  may  be  fixed  with  a  specified 
beamwidth  (ground  fixed),  or  continually  tracking  the  vehicle  with  a  specified 
beamwidth  (ground  tracking),  or  it  may  be  isotropic.  There  are  nine  possible 
comi)inations  of  these  antennas,  but  the  vehicle-tracking  isotropic-ground  case 
is  not  computed  here.  However,  the  performance  of  this  combination  can  be 
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obtained  as  a  special  case  of  the  vehicle-tracking  fixed-ground  case  by  choosing 
the  proper  beamwidth  for  the  ground  antenna.  Space  losses  are  also  printed  out 
separately. 

In  subroutine  RUNG  a  numerical  integration  of  the  true  anomaly  equation  is 
performed  by  a  fourth-order  Runge-Kutta  method. 

There  are  no  READ  or  PRINT  statements  in  either  DEBE  or  RUNG. 

Subroutines  TRAC  and  RUNG  together  solve  the  equation  of  motion  of  the 
vehicle  and  generate  the  ephemeris  data  for  the  vehicle.  The  use  of  CALC  and 
DEBE  is  optional  depending  on  the  application  under  investigation.  Thus,  the 
over-all  program  is  made  up  of  engineering  building  blocks  which  can  be  put 
together  to  solve  a  particular  problem  with  economy  of  computer  time  and  a 
minimum  of  extraneous  information  in  the  output.  A  further  advantage  of  this 
approach  is  that  refinements  and  additions  can  be  incorporated  into  the  individual 
subroutines  without  the  necessity  of  recompiling  the  whole  program. 

INPUT  DATA  CARDS 

The  proper  sequencing  of  cards  for  the  running  of  multiple  cases  is 
described  under  the  heading  Stacking  Cases.  The  program  is  set  up  to  determine 
the  orbital  coverage  capabilities  of  a  given  arbitrary  set  of  global  tracking  stations 
for  n-different  orbit  configurations.  The  number  of  input  data  cards  in  any  given 
run  will  depend  on  both  the  number  of  stations  and  the  number  of  orbits  investi¬ 
gated.  However,  there  are  only  four  basic  types  of  input  data  cards. 

Card  #1  contains  the  burnout  data  which  establishes  the  orbit.  Read-in 
occurs  in  program  MAIN.  There  will  be  one  such  card  for  each  orbit  configuration 
investigated.  The  second  to  n-sets  of  burnout  data  cards  are  arranged  sequen¬ 
tially  at  the  end  of  data  input  deck. 
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Card  #2  contains  the  control  constants  which  are  used  by  the  program. 

These  constants  control  the  integration  interval,  printout  intervals  and  increments, 
and  the  subroutine  and  the  time-stop  options.  The  printout  interval  Is  the  sub¬ 
division  of  time  in  which  the  data  is  printed  out  within  the  larger  printout 
increment.  The  purpose  of  the  printout  increment  option  is  to  allow  for  a  "fine¬ 
grained"  printout  during  portions  of  the  trajectory  that  are  changing  rapidly 
(or  are  of  more  interest  than  some  other  part  of  the  trajectory),  and  conversely, 
to  allow  a  "coarse-grained"  printout  during  slowly  varying  parts  of  the  trajectory. 
The  unused  control  constants  were  inserted  for  use  in  furture  program  modifica¬ 
tions.  Only  one  such  card  will  appear  in  any  given  machine  run.  Read-in  occurs 
in  subroutine  TRAC. 

Card  #3  contains  ground  and  vehicle  sensor  data.  Only  one  such  card  will 
appear  in  any  given  machine  run.  Read-in  occurs  in  subroutine  CALC. 

Card  #4  contains  the  geographic  location,  minimum  antenna  elevation  angle, 
and  beamwidth  data  of  the  tracking  stations.  There  will  be  one  card  for  each 
station,  and  these  cards  will  be  arranged  sequentially  when  more  than  one  statio! 
is  being  investigated.  Read-in  occurs  in  subroutine  CALC. 

The  card  formats  and  Fortran  nomenclature  are  summarized  in  Tables  1, 

2,  3,  and  4. 

STACKING  CASES 

For  runs  with  one  orbit  configuration  and  multiple  ground  tracking  stations, 
the  proper  data  card  stacking  is: 

Card  #1,  burnout  conditions 

Card  #2,  control  constants 

Card  #3,  vehicle  and  ground  sensor  data 
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Table  1 


Burnout  Data 


CARD  #1  (One  Card) 


Fortran 

Variable 

Card 

Column 

Location 

Fortran 

Numeric 

Format 

Definition 

HBO 

1-10 

F10.2 

Burnout  altitude  -  measured  positive 
UP  from  surface  of  earth  (NM) 

VBO 

11-20 

F10.2 

Burnout  velocity  (FT/SEC) 

GBO 

21-30 

F10.2 

Burnout  flight  path  angle  -  measured 
positive  up  from  local  horizontal 
plane  (DEG) 

BOHEAD 

31-40 

F10.2 

Burnout  heading  -  measured  positive 
clockwise  from  north  (DEG) 

BOLAT 

41-50 

F10.2 

Burnout  latitude  -  measured  positive 
North  from  the  Equator  (DEG) 

BOLONI 

51-60 

F10.2 

Burnout  longitude  -  measured  positive 
East  from  Greenwich  (DEG) 

BIG 

61-70 

F10.2 

Burnout  spin  axis  attitude  -  measured 
positive  up  from  local  horizontal 
plane  (DEG) 

BIH 

71-80 

F10.2 

Burnout  spin  axis  heading  -  measured 
positive  clockwise  from  North  (DEG) 
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Table  2 


Program  Control  Constants 


Card  #2  (One  Card) 


Fortran 

Variable 

Card 

Column 

Location 

Fortran 

Numeric 

Format 

Definition 

NUM  (1) 

1-5 

15 

Number  of  ground  sensors  (maximum 
of  150  sensors) 

NUM  (2) 

6-10 

15 

Number  of  integration  intervals  in  the  | 
initial  and  final  printout  increment 

NUM  (3) 

11-15 

15 

End  of  first  printout  increment  in  ; 

seconds 

....  ... 

NUM  (4) 

16-20 

15 

Number  of  integration  interval  in 
intermediate  printout  increment  ' 

NUM  (5)** 

21-25 

15 

Final  time  in  seconds  or  as  modified 
by  NUM  (16) 

NUM  (6) 

26-30 

15 

Control  constant  to  select  time  stop 
option 

^  0  bypass  option  -  use  increments 
>  0  select  option  . 

NUM  (7) 

31-35 

15 

Control  constant  to  select  CALC- 
subroutine  option  ' 

£  0  bypass  option  -  use 
>  0  select  option 

NUM  (8) 

36-40 

15 

Integration  interval  in  seconds 

NUM  (9) 

41-45 

15 

Control  constant  to  select  DEBE  - 
subroutine  option 
>  0  bypass  option 
^  0  select  option 

NUM  (16)* 

76-80 

15 

Constant  which  multiplies  NUM  (5) 
to  extent  final  time;  i.  e.  ,  when 

NUM  (16)  =  1  ,  NUM  (5)  is  in  sec. 

NUM  (16)  =  60  ,  NUM  (5)  is  in  min. 

NUM  (16)  =  3600,  NUM  (5)  is  in  hrs, 

etc. 

*NOTE:  NUM  (10)  through  NUM  (15)  are  unassigned. 

**NOTE:  In  the  program  printout  of  the  program  control  constants  NUM  (5) 
reverts  to  the  total  program  time  in  seconds. 
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Table  3 


Sensor  Data  -  Ground  and  Vehicle 


Card  #3  (One  Card) 


Fortran 

Variable 

Card 

Column 

Location 

Fortran 

Numeric 

Format 

Definition 

FREQ 

1-10 

F10.2 

Antenna  operating  frequency  in 
megacycles 

BWl 

11-20 

F10.2 

Beamwidth  of  spaceborne  antenna 
measured  between  half-power  points  in 
degrees 
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Table  4 


Station  Data 


Card  #4  (NUM  (1)  Cards*) 


Fortran 

Variable 

Card 

Column 

Location 

Fortran 

Format 

Definition 

9-24 

2A8 

These  16  characters  reserved  for  the 

station  name 

ALAT  (I) 

25-32 

F8.  2 

Latitude  of  the  Ith  ground  sensor 
measured  positive  north  of  equator 
(Deg.) 

ALONG  (I) 

33-40 

F8.2 

Longitude  of  the  Ith  ground  sensor  - 
measured  positive  east  of  Greenwich 
meridian  (Deg. ) 

ANEL  (I) 

41-47 

F8.2 

Elevation  of  the  Ith  ground  sensor,  if 
fixed  -  measured  positive  above 
tangent  plane  (Deg. ) 

ANAZ  (I) 

48-56 

F8.2 

Azimuth  of  the  Ith  ground  sensor,  if 
fixed  -  measured  positive  clockwise 
from  north  (Deg. ) 

ALTS  (I) 

57-64 

F8.2 

Altitude  of  station  above  seal  level 
(feet) 

BW  (I) 

65-72 

F8.2 

Beamwidth  of  the  Ith  ground  sensor  - 
measured  between  half-power  points 
(Deg. ) 

ELEM  (I) 

73-80 

F8.  2 

Minimum  elevation  angle  of  tracking 
antenna  (Deg. ) 

*  One  of  these  cards  is  required  for  each  station.  The  Fortran  variable  NU^l(l) 
mu.st  coincide  with  the  number  of  station  cards.  The  number  of  stations  can 
vary  from  1  to  150. 


13 


Card  #4,  station  1 


t  t 

Card  #4,  station  n 


For  runs  with  more  than  one  orbit  configuration  and  ground  tracking  station, 
the  proper  data  card  stacking  is 

Card  #1,  first  burnout  conditions 

Card  #2 ,  control  constants 

Card  #3,  vehicle  and  antenna  sensor  data 


Card  #4,  station  1 


Card  #4,  station  2 


I  I 

;  I 

Card  #4,  station  n 


Card#l,  second  burnout  condition 


I  I 

I  I 

I  I 


t  t 

Card  #1,  mth  burnout  conditions 


OUTPUTS 


In  the  most  general  case,  when  all  functional  options  are  selected,  the 
program  delivers  the  following  five  types  of  output  data;  burnout  conditions, 
orbital  parameters  at  burnout,  station  data,  vehicle  and  sensor  parameters,  and 

coverage  data.  *  The  coverage  data  includes  the  vehicle  ephemeris  data  and  the 
performance  data  for  each  station  that  is  in  the  geometrical  line  of  sight  of  the 


*  For  convenience  in  checking  the  options  that  have  been  selected,  the  program 
control  constants  as  defined  in  Table  2  appear  at  the  beginning  of  the  printout. 
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vehicle  as  a  function  of  time.  If  the  CALC  option  is  not  selected,  only  the 
vehicle  burnout  and  ephemeris  data  will  be  printed.  The  definition  of  each 
variable  in  the  five  groups  as  it  appears  on  the  printed  output  is  described  in 
Tables  5  through'9.  A  typical  printout  is  shown  in  Table  10. 


Table  5 

Burnout  Conditions 


Output  #1  * 


Variable 

Units 

Definition 

ALTITUDE 

n.  m. 

Burnout  altitude  from  the  surface 
of  the  earth 

VELOCITY 

ft/ sec 

Inertial  burnout  velocity 

GAMMA 

deg. 

Inertial  burnout  flight  path  angle 

HEADING 

deg. 

Inertial  heading  at  burnout 

LATITUDE 

deg. 

Vehicle  latitude  at  burnout 

LONGITUDE 

deg. 

Vehicle  longitude  at  burnout 

SPIN  AXIS  ATTITUDE 

deg. 

Antenna  axis  attitude  at  burnout 

SPIN  AXIS  HEADING 

deg. 

Inertial  heading  of  the  antenna  axis 
at  burnout 

*  These  outputs  correspond  exactly  to  the  burnout  data  read  in  on  card  #1 
(see  Table  II). 
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Table  6 


Orbital  Parameters  At  Burnout 


Output  #2 


Variable 

Units 

Definition 

ECCENTRICITY 

none 

Eccentricity  of  the  orbit 

INCLINATION 

deg. 

Inclination  of  the  orbital  plane 
<90^  prograde,  >  90^  retrograde 

SEMI  LAT  REC 

n.  m. 

Semi-latus  rectum 

TRUE  ANOMALY 

deg. 

Angle  between  perigee  and  the 
vehicle  in  the  orbital  plane  at 
burnout 

OMEGA 

deg. 

Angle  between  ascending  node  and 
perigee  in  the  orbital  plane  at 
burnout 

ASNODE 

deg. 

Angle  between  the  ascending  node 
and  Greenwich  (+  East) 
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Table  7 


Station  Data 


Output  #3* 


Variable 

Units 

Definition 

NUMBER 

none 

Station  number 

NAME 

none 

Station  name  limited  to  16  characters 

LAT 

deg. 

Station  latitude  (+  north  of  the 

Equator) 

LONG 

deg. 

Station  longitude  (+  east  of  the 
Greenwich  meridian) 

ELEV 

deg. 

Elevation  angle  of  the  fixed  ground 
sensor  axis  above  the  tangent  plane 

AZIMUTH 

deg. 

1 

1 

Azimuth  of  the  fixed  ground  sensor  | 
axis  (+  clockwise  from  north) 

ALTITUDE 

feet 

Altitude  of  the  station  above  sea 
level 

BEAMWIDTH 

deg. 

Beamwidth  of  the  ground  sensor 
measured  between  half  power  p>oints 

MIN  ELEV  ANGLE 

deg. 

Minimum  elevation  angle  of  ti’acking 
antenna 

*  These  outputs  correspond  to  the  station  data  read  in  on  card  #4  (one  for  each 
station)  (see  Table  4). 
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Table  8 


Vehicle  And  Sensor  Parameters 


Output  #4 


Variable 

Units 

Definition 

BEAMWIDTH 

deg. 

Vehicle  antenna  beamwidth  measured 
between  the  half-power  points 

TRACKING  GAIN 

db 

Vehicle  on  axis  antenna  gain 

OPERATING  FREQUENCY 

megacycles 

Operating  frequency  of  the  sensors 

NOTE:  The  definition  of  the  losses  for  the  various  vehicle-ground  antenna 

configurations  appears  following  the  above  output  for  easy  reference 


Table  9 
Coverage  Data 

Output  #5  Part  I  -  Vehicle  Position 


Variable 

Units 

Definition 

TIME 

min. 

Instantaneous  time  in  orbit 

TIME 

hrs. 

Instantaneous  time  in  orbit 

LATITUDE 

deg. 

Vehicle  latitude  in  orbit 

ALTITUDE 

min. 

Vehicle  altitude  above  the  surface  of 
the  earth 

TRUE  ANOMALY 

deg. 

True  anomaly 

OMEGA 

deg. 

Angle  between  perigee  and  ascending 
node  in  the  orbital  plane 

ASNODE 

deg. 

Angle  between  the  ascending  node 
and  the  Greenwich  meridian 
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Table  9  (Cont' d) 


Part  II  -  Stations  in  Contact 


Variable 

Units 

Definition 

ELEVATION* 

deg. 

Elevation  angle  of  the  line  of  sight 
from  the  station  to  the  vehicle 
{+  up  from  the  tangent  plane) 

AZIMUTH* 

deg. 

Azimuth  angle  of  the  line  of  sight 
between  the  station  and  the  vehicle 

RANGE* 

n.  m. 

Line  of  sight  distance  between  the 
station  and  the  vehicle 

ANT  ANGLE* 

SPIN  STAB 

deg. 

Angle  between  spin-stabilized  vehicle 
antenna  and  the  line  of  sight  between 
the  station  and  the  vehicle 

ANT  ANGLE* 

SPIN  STAB 

deg. 

Angle  between  the  earth-centered 
vehicle  antenna  and  the  line  of  sight 
between  the  station  and  the  vehicle 

GRD  ANT* 

ANGLE 

deg. 

Angle  between  the  fixed  ground 
antenna  axis  and  the  line  of  sight 
between  the  station  and  the  vehicle 

VEHICLE  ANT 
GAINS 

db 

Gains  of  the  space-borne  antennas 

SPIN  STAB 

db 

Spin-stabilized  antenna 

VERT-STAB 

db 

Earth-centered  antenna  gain 

GRD  ANT  GAIN 

db 

Gain  of  the  fixed  ground  antenna 

SPACE  LOSS 

db 

Total  space  loss 

DBl  through  DBS 

db 

Total  losses  for  the  various  vehicle 
ground  antenna  configurations. 
Definitions  appear  on  each  printout 
(see  DEBE  listing  in  Appendix  E) 

NOTE:  The  rate  of  change  of  each  of  the  parameters  marked  with  an  asterisk  is 
printed  out  directly  below  the  value  of  each  of  the  parameters. 
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POWERED  FLIGHT  OPTION 


The  complete  engineering  simulation  of  the  space-ground  environment 
covers  the  entire  mission  of  the  space  vehicle  from  launch  through  the  orbital 
motion.  In  many  situations,  however,  the  detailed  solution  from  launch  is  not 
required.  Program  MAIN  permits  the  convenient  division  of  the  simulation  at 
the  point  of  injection. 

The  over-all  program  can  be  operated  in  two  basic  modes: 

In  the  first,  or  normal  mode,  the  complete  solution  from  launch  through 
orbital  motion  is  obtained,  and  program  MAIN  is  not  used.  Burnout  conditions 
are  calculated  directly  in  the  powered-flight-trajectory  program  described  in 
Reference  [l]  and  passed  on  to  the  proper  subroutines  which  compute  orbital 
position  and  station  coverage.  * 

In  the  second  mode,  the  burnout  conditions  are  not  calculated  but  inserted 
into  the  program  on  input  data  card  #1.  Thus,  the  launch-trajectory  portion  of 
the  simulation  is  not  used.  Program  MAIN  then  functions  as  the  executive 
program  for  the  reading  and  routing  of  these  parameters  into  the  subroutines 
computing  orbital  position  and  station  coverage.  The  program  also  provides  for 
the  proper  end-of-program  statements  for  the  logical  Fortran  IV  termination  of 
the  program. 

SPECIAL  ROUTINES 

There  are  only  two  special  subroutines  required  to  run  the  program  on  the 
IBM  7030  computer.  These  are  the  inverse  sine  and  cosine  functions  of  an  angle. 
These  routines  were  prepared  for  use  with  this  simulation,  and  their  Fortran  IV 
lists  are  included  in  Appendix  VI. 

*  The  only  modification  required  in  the  Powered  Flight  Trajectory  Program  as 
described  in  Reference  [l]  is  the  addition  of  a  CALL  TRAC  (...)  statement 
following  the  test  for  last- stage  burnout. 
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SECTION  III 


ORBITAL  SIMULATION  -  SUBROUTINE  TRAC 


Subroutine  TRAC,  called  by  program  MAIN  or  the  powered-flight-trajectory 
program,  has  two  primary  functions.  First,  it  calculates  the  time-varying  orbital 
position  of  the  vehicle  relative  to  geocentric  axes  fixed  in  a  spherical  rotating 
earth.  Second,  it  performs  the  basic  control  functions  of  CALC  subroutine 
selection  and  of  integration  interval  and  printout  interval  management.  A  func¬ 
tional  flow  diagram  of  TRAC  is  given  in  Fig.  3  and  its  Fortran  IV  listing  is  given 
in  Appendix  II. 

CALCULATION  OF  FIXED  ORBITAL  PARAMETERS 

Orbital  position  determination  as  a  function  of  time  starts  with  the  calcu¬ 
lation  of  the  fixed  orbit  geometry  and  orbit  position  relative  to  the  earth  from 
the  burnout  conditions.  The  engineering  simulation  described  in  this  report 
assumes  a  spherical-earth  model  and  neglects  perturbation  effects,  such  as  nodal 
precession  and  rotation  of  the  line  of  perigee  in  the  orbital  plane.  Consequently, 
once  the  burnout  conditions  are  specified,  namely,  altitude  (HBO),  velocity  (VBO), 
flight-path  angle  (BIG) ,  heading  (BOHEAD) ,  latitude  (BOLAT) ,  and  longitude 
(BOLONG),  the  position  of  the  orbital  plane  and  the  orbit  of  the  vehicle  in  that 
plane  are  completely  established  and  fixed  in  time.  This  geometry  is  completely 
described  by  the  following  fixed  orbital  parameters  which  are  computed  in  sub¬ 
routine  TRAC. 

(1)  PLAIN  -  the  orbital  plane  inclination  relative  to  earth. 

(2)  VEE  -  initial  true  anomaly  or  the  angle  between  the  line  of  perigee  and 
the  position  of  the  vehicle  in  orbit  in  the  orbital  plane  at  injection. 
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ENTER  SUBROUTINE 


PROGRAM 


Fig.  3  Functional  Flow  Diagram — Subroutine  TRAC 
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(3)  ECC  -  eccentricity  of  the  orbit. 

(4)  P  -  semi-latus  rectum. 

(5)  OMEGA  -  perigee  position  relative  to  the  ascending  node  in  the  plane 
of  the  orbit. 

(6)  ASNODE  -  longitude  of  ascending  node  relative  to  an  earth  reference 
meridian. 

These  quantities  are  shown  in  Fig.  4.  The  mathematical  expressions  for  each 
of  these  quantities  are  derived  below. 

Inclination  Angle  (PLAIN) 

The  inclination  angle  can  be  readily  determined  from  the  spherical  triangle 
shown  in  Fig.  4. 

From  the  cosine  law 


cos  PLAIN 


cos  BC  -  cos  AC  cos  AB 
sin  AC  sin  AB 


(1) 


cos  AB 


cos  BC  cos  AC  -  sin  BC  sin  AC  cos  C 


cos  DC  cos  AC  (since  C  =  90°) 


(2) 


From  the  sine  law 


sin  AB  _  sin  AC 
sin  C  sin  B 


or  sin  AB  =  sin  AC/ sin  B  (since  C  =  90°) 


(3) 
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B 


SPHERICAL  TRIANGLE  ABC 


Fig.  4 


Orbital  Geometry 


Substituting  Eqs.  (2)  and  (3)  in  Eq.  (1)  gives 


cos  (PLAIN) 


cos  BC  -  cos  BC  cos  AC  cos  AC 
sin^  AC/ sin  B 


2 

=  cos  BC  (1  -  cos  AC)  sin  B 
2 

sin  AC 

=  cos  BC  sin  B 


(4) 


where  BC  is  the  burnout  latitude 

B  is  the  burnout  heading  angle 

In  Fortran  notation 

PLAIN  =  ACOS  (CLABO  *  SHEBO  (5) 

where 

*  =  multiplication 

CLABO  =  cosine  of  burnout  latitude 

SHEBO  =  sine  of  the  burnout  heading  angle.  If  SHEBO  is  positive,  i.e.  ,  the 
burnout  heading  between  0  and  180  degrees,  PLAIN  will  be  between 
0  and  90  degrees,  or  a  prograde  orbit.  If  SHEBO  is  negative,  i.e.  , 
the  burnout  heading  between  180  and  360  degrees,  then  PLAIN  will  be 
between  90  and  180  degrees,  or  a  retrograde  orbit.  The  term  CLABO 
will  always  be  positive. 

ACOS  =  arc  cosine 
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True  Anomaly  at  Burnout  (VEE) 


The  angle  between  the  perigee  and  the  point  in  the  orbit  at  burnout  is  the 
initial  true  anomaly  of  the  ellipse  and  is  given  by  the  following  relation. ' 


VEE  =  true  anomaly  =  tan 


(HBO)  (VBO)^ 

GMU 


sin  (GBO)  cos  (GBO) 


(ALT)  (VBO)^ 

GMU 


2 

cos 


(GBO)  -1 


(6) 


where  GMU  is  the  gravitational  mass  unit  or  constant  K  in  Reference  [2] .  The 
magnitude  of  the  true  anomaly  depends  on  the  signs  of  the  numerator  and  denom¬ 
inator,  both  of  which  can  be  plus  or  minus  independently  of  the  other.  In  the 
numerator,  the  burnout  flight  path  angle  GBO  can  be  plus  or  minus  relative  to  the 
local  horizon.  The  term  reflecting  this  sign  is  sin  (GBO)  which  can  be  either 
plus  or  minus.  In  the  denominator,  the  first  term  can  be  greater  or  less  than 
one.  The  sign  is  primarily  determined  by  the  burnout  velocity. 

The  simulation  logic  must  be  arranged  to  provide  the  proper  value  of  the 
true  anomaly  for  different  combinations  of  signs  of  the  numerator  and  demon.- 
inator  at  burnout.  A  geometric  interpretation  of  these  different  combinations 
is  given  in  Fig.  5.  Detailed  Fortran  logic  of  this  part  of  the  progp'am  is  given 
on  lines  T38  through  T41  in  Appendix  II. 


Eccentricity  (ECC) 


The  eccentricity  of  the  ellipse  is  given  by  the  equation: 


I2] 


ECC  = 


(HBO)  (VBO)^  -  1 

cos^  (GBO)  +  sin^  (GBO) 

GMU 

(7) 
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.  (HBOXVBO)* 
A  « - 

G  MU 


CASE  I  Acoi*  GBO  -  I  >0 
SlnGB0>0 


Bo  *  tan-' 


(+} 

(+) 


VEE  >  Bo 


CASE  2  A  coi*  GBO  -  I  >0 


Sin  GB0<0 


»<.  ■  tan-i  jij 

VEE  =  360- I ©o I 


CASE  3  A  CO**  GBO- I  <0  Bo  =  TAN*' 

Sin  GB0>0  a-TAN"' 

VEE  *  ISO  -  |9ol 

*  ISO  +  ©0 


CASE  4  Aco**GBO-l<0  Bo  *  TAN"' 

Sin  GB0<0  =  TAN"' 

VEE  =  ISO  +  I  ©o  I 
=  ISO  +  Bo 


Fig.  5  Quadrant  Determination  of  True  Anomaly  at  Burnout 
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where 


HBO  =  altitude  in  nautical  miles  from  earth  center  to  burnout  altitude 
VBO  =  burnout  velocity  in  knots 

2  3  2 

GMU  =  gravitational  mass  unit  =  gR^  =  62 , 750.  21  n.  m.  per  sec. 

GBO  =  burnout  flight  path  angle 
Semi-latus  Rectum  (P) 

The  distance  from  the  focal  point  of  the  ellipse  to  a  point  on  the  elliptical 
orbit  where  the  true  anomaly  is  equal  to  90  degrees  is  called  the  semi-latus 
rectum  or  the  parameter  of  the  orbit  (the  line  FB  in  Fig.  4. )  It  is  given  by  the 
polar  equation  for  a  conic.  ^ 

P  =  HBO  (1  +  ECC  *  cos  VEE)  (8) 

where  these  quantities  have  all  been  defined  previously. 

Argument  of  Perigee  at  Burnout  Relative  to  the  Ascending  Node  (OMEGA) 

The  argument  of  perigee,  OMEGA  in  Fortran  format,  is  generated  from 
the  true  anomaly  (VEE)  and  the  arc  AB  of  the  spherical  triangle  ABC  of  Fig.  4. 
The  relations  which  exist  between  OMEGA,  true  anomaly,  and  arc  AB  depend 
upon  the  burnout  heading  and  latitude. 

In  Fig.  6,  the  possible  permutations  of  these  two  burnout  parameters  and 
the  corresponding  expressions  for  the  argument  of  perigee  OMEGA  are  given. 

For  any  given  burnout  heading  quadrant,  it  should  be  noted  that  changing  from 
a  northern  to  southern  hemisphere  burnout  location  merely  changes  the  sign  of 
the  arc  AB  in  the  expression  for  OMEGA  (AP).  Recognition  of  this  characteristic 
simplifies  the  Fortran  program  logic. 


34 


(q)NORTHERN  hemisphere 
(Sin(LAT)  >0  :Cot(LAT)  >0) 


CASE  I  CASE  2 


CASE  3  CASE  4 


Co5(HBO|>0:  Sln(HB0)>0  Cos(HBO|<0:S(n  (HBO)>0 


A 


ORBITAL  PLANE 
[■BURNOUTl^ 

'®°^LheadingJ-^^° 

Cot(HBOI<0:  Sin  (HBO)  <  0 


EQUATORIAL  PLANE 


ORBITAL  PLANE 
[BURNOUT] < 

^  [heading 

Cos(HBO)  >0:  Sin  (HBO)  <0 


APs-PB-l-AB 

AP:  IT-PB-A'B 

AP  :  it-PB-a'  B 

APs-PB-l-AB 

XA=  XC- AC 

XA^ir  +XC+A'C 

XA  :  Tr+XC±A'C 

XA=XC  +  AC 

WHERE  : 


PB  =  true  ANOMOLY  (VEE) 

AP  =  ARGUMENT  OF  PERIGEE  RELATIVE  TO  ASCENDING  NODE  (OMEGA) 

XA  =  LONGITUDE  OF  ASCENDING  NODE  (ASNOOE) 

XC  -  LONGITUDE  OF  BURNOUT  (BOLONG) 

AB  ARGUMENT  OF  BURNOUT  RELATIVE  TO  ASCENDING  NODE  IN  ORBIT  PLANE 
AC  =  ARGUMENT  OF  BURNOUT  RELATIVE  TO  ASCENDING  NODE  IN  EQUATORIAL  PLANE 


Fig.  6  Quadrant  Resolution  of  Argument  of  Perigee  and  Ascending  Node 
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(b) SOUTHERN  HEMISPHERE 
(Sin(LAT)<0  :Cot(LAT)>0) 


EQUATORIAL  PLANE  EQUATORIAL  PLANE 


ORBITAL  PLANE  ORBITAL  PLANE 


ORBITAL  PLANE  ORBITAL  PLANE 


[BURNOUT 

Lheading 


]<- 


<90 


90 


BURNOUT 

HEADING 


<  180 


I80< 


BURNOUT 

.HEADING 


<270 


270  < 


BURNOUT 

.HEADING 


<360 


AP=-PB-AB 

AP=7r-PB+A'B 

AP  =  Tr-PB+A'B 

AP  =  -PB^-AB 

XA=  XC  +  AC 

XA=7r+XC-/<C 

XA  s  IT  +XC^-A'C 

XA*XC  -  AC 

Fig.  6  (Continued)  Quadrant  Resolution  of  Argrument  of  Perigee  and 
Ascending  Node 
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The  trigonometrical  relations  for  the  evaluation  of  AB  follows  directly 
from  the  sine  law.  Referring  to  Fig.  4,  sin  C  =  sin  90°  =  1 


AB 


ARC  sin 


sin  BC 
sin  (PLAIN) 


ARC  sin 


sin  (latitude) 

sin  (inclination) 


(9) 


or  in  Fortran  notation 

AB  =  ASIN  (SLAIN)*  (10) 

where  SLAIN  =  SLABO/SINL 

The  Fortran  listing  of  the  branching  logic  for  this  portion  of  the  program  is 
contained  on  lines  T54  through  T61. 1  in  Appendix  II. 

Longitude  of  the  Ascending  Node  (ASNODE) 

The  longitude  of  the  ascending  node  is  called  ASNODE  in  the  Fortran 
program.  Its  value  also  depends  on  the  burnout  heading  and  latitude.  Prograde 
and  retrograde  launchings  will  generate  ascending  nodes  differing  by  180  degrees 
due  to  the  definition  of  the  ascending  node.  The  relations  for  the  longitude  of 
the  ascending  node  are  given  in  Fig.  6  for  the  possible  burnout  latitude  and 
heading  permutations.  Similar  to  the  argument  of  perigee,  changing  from  a 
northern  to  a  southern  burnout  location  at  a  given  heading  merely  changes  the 
sign  of  the  arc  AC  in  the  expression  for  the  ASNODE  (arc  XA,  Fig.  6). 


*  The  ASIN  subroutine  treats  a  negative  angle  as  minus  the  arc  sine  of  the 
magnitude  of  the  angle.  This  is  consistent  with  the  branching  logic,  since 
AB  will  be  positive  for  the  Northern  Hemisphere  and  negative  for  the 
Southern  Hemisphere. 
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The  trigonometrical  expression  for  arc  AC  in  terms  of  known  orbital 
data  can  again  be  derived  directly  from  the  basic  spherical  triangle  in  Fig.  4 
and  the  sine  law. 

sin  C  =  1  .  .  sin  AC  =  sin  AB  sin  B  (11) 


Substituting  for  sin  AB  from  Eq.  (9) 


sin  AC  = 


sin  BC  sin  B 
sin  (PLAIN) 


but  B  is  the  burnout  heading;  thus,  in  Fortran  notation, 


AC  =  ASIN  (SLABO  *  SHEBO/SINL) 
or  AC  =  ASIN  (SLAIN  *  SHEBO) 


(12) 


(13) 


where  SHEBO  is  the  sine  of  the  burnout  heading  angle. 

The  corresponding  Fortran  listing  of  the  branching  logic  is  contained  on 
lines  T54  through  T61. 1  in  Appendix  II. 

CALCULATION  OF  TIME -VARYING  ORBITAL  PARAMETERS 

The  fixed  parameters  given  above  establish  the  initial  position  of  the 
vehicle  relative  to  the  earth.  From  these  six  initial  conditions  and  the  differ¬ 
ential  equation  of  motion  for  the  variation  of  the  true  anomaly  with  time,  the 
position  of  the  vehicle  relative  to  the  moving  earth  can  be  determined.  The 
variation  of  the  true  anomaly  is  obtained  by  a  numerical  integration  of  the 
equation  of  motion  of  the  orbiting  vehicle.  This  calculation  is  performed  in 
subroutine  RUNG  whose  calling  is  controlled  by  subroutine  TRAC.  Subroutine 
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RUNG  then  outputs  the  true  anomaly  to  subroutine  TRAC  where  the  calculation 
of  the  following  three  parameters  defining  the  vehicle  position  relative  to  a 
moving  earth  are  made. 

RLAT  =  Latitude  of  vehicle  on  moving  earth. 

RLONG  =  Longitude  of  vehicle  on  moving  earth. 

ALT  =  Altitude  of  vehicle  above  spherical  earth. 

The  calculation  of  the  time -varying  orbital  parameters  begins  with  statement 
number  700  of  the  TRAC  listing  in  Appendix  II. 

Latitude  of  Vehicle  on  Moving  Earth  (RLAT) 

The  rotation  of  a  spherical  earth  has  no  net  effect  on  the  latitude  of  a 
vehicle  relative  to  the  earth.  An  expression  relating  the  latitude  at  time  t.  and 
the  parameters  previously  given  can  be  readily  obtained  from  the  spherical 
triangle  ABC  in  Fig.  4. 

From  the  law  of  sines 

sin  (RLAT)  =  sin  BC  =  sin  (PLAIN)  sin  AB  (14) 

AB  =  VEE  +  OMEGA  (15) 

Substituting  into  Eq.  (14)  and  expanding  the  sine 

sin  (RLAT)  =  sin  (PLAIN)  [sin  (VEE)  cos  (OMEGA)  +  cos  (VEE)  sin  OMEGA] 

(16) 
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In  Fortran  listing  the  following  notation  is  used. 


sin  (RLAT)  =  SLAR;  sin  (VEE)  =  SVEE;  sin  (OMEGA)  = 

sin  (PLAIN)  =  SINL;  cos  (VEE)  =  CVEE;  cos  (OMEGA)  =  COM 

Cg  =  SOM  *  SINL 

C,  =  COM  *  SINL 
4 

where  VEE  is  now  the  true  anomaly  at  time  t. 


1 


SLAR  =  C  *  SVEE  +  C  *  CVEE 

4  O 


(17) 


RLAT  =  ASIN  (SLAR)  *  CONVI 


(18) 


where  CONVI  is  a  conversion  factor  from  radians  to  degrees.  The  geometric 
relationships  for  OMEGA  and  VEE  were  given  in  Fig.  6. 

Longitude  of  Vehicle  on  Moving  Earth  (RLONG) 

The  earth' s  rotation  will  affect  the  longitude  of  the  vehicle  to  the  earth, 
since  the  reference  from  which  the  longitude  is  measured  is  fixed  in  the  earth. 
Let  WE  be  the  angular  velocity  of  the  earth;  then  at  time  t.,  the  longitude  is 
given  by  the  expression: 


RLONG  (t.)  =  CONGR  =  ASNODE  -  WE  (t.)  ±  AC 


(19) 


where  from  the  cosine  law  AC  is  given  by  the  following  expression: 


cos  (OMEGA  +  VEE  (t.)) 


(20) 


cos  (LAT) 
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Thus,  in  terms  of  previously  defined  Fortran  notation 


CONGR  =  ASNODE  ±  ACOS  (COM  *  CVEE  -  SOM  *  SVEE)/(CLAR)-WE  *  TIME 

(21) 

and  RLONG  =  CONGR  *  CONVI  (22) 

A  plus  sign  with  the  AC  term  in  the  expression  CONGR  signifies  a  prograde 
orbit;  a  minus  sign  signifies  a  retrograde  orbit.  The  burnout  heading  angle 
which  determines  whether  the  motion  is  prograde  or  retrograde  can  be  used  to 
control  the  branch  logic.  This  logic  is  contained  on  lines  T82  through  T87  in 
Appendix  II. 

The  quadrant  of  the  term  cos  (OMEGA  +  VEE  (t. ))  in  the  above  expression 
for  AC  is  lost  when  this  angle  exceeds  n.  That  is,  the  particular  inverse 
cosine  subroutine  used  in  the  program  only  gives  angle  values  between  0  and  n. 
This  programing  difficulty  was  circumvented  by  using  the  sine  of  the  latitude 
(SLAR)  to  control  the  branching  to  either  the  plus  or  minus  form  of  the  above 
expression  for  CONGR.  By  definition  of  the  ascending  node,  the  angle  (OMEGA 
+  VEE)  is  between  n  and  2n  in  the  Southern  Hemisphere.  The  nature  of  these 
relationships  are  illustrated  in  Fig.  7. 

Orbital  Altitude  (ALT) 

For  a  spherical -earth  model,  the  only  time-varying  parameter  in  the 
expression  for  orbital  altitude  is  the  true  anomaly,  or 

ALT  =  P/(1.0  +  cos  (VEE(tp)  (23) 

where 

P  =  the  semi-latus  rectum 
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C, 


A  =  OMEGA 

WHEN  n-r  <  +  VEE  <  (n  + I  )ir  :  CASE  C, 

LATITUDE  >  0 

CONOR*  AN-Wi TIME  +  CO«-'  (co«  180-0  ) 
«AN-Wi  TIME  +(180-0) 

WHEN  (n +I)T  <  +  VEE  <  (n  +  2  IT  :  CASE  C2 

latitude  <  0 

=  AN-We  time  -cot*'  (cot  180  +  0) 
»  AN  -We  time  -H 


Fig.  7  Resolution  of  Longitude  Ambiguity 

Note  that  this  is  the  same  as  Eq.  (8)  except  that  VEE  is  the  true  anomaly  at 

time  t.  rather  than  at  burnout. 

1 

DIRECTION  COSINES  OF  THE  SPIN-STABILIZED  ANTENNA  (XI,  YI,  ZI) 

The  calculation  of  the  direction  cosines  of  the  spin-stabilized  antenna  axis 
relative  to  the  X,  Y,  Z  earth-centered  coordinates  are  made  in  subprogram 
TRAC.  This  represents  a  slight  deviation  from  the  basic  approach  of  having 
each  subroutine  perform  a  self-contained  engineering  function.  These  direction 
cosines  are  passed  to  subroutine  CALC  for  computing  the  geometrical  relation¬ 
ships  between  the  orbital  vehicle  antenna  axis  and  the  ground  sensors.  These 
direction  cosines  involve  the  same  functions  of  the  burnout  conditions  as  the 
orbital  parameters  at  burnout;  therefore,  by  making  the  calculation  in  TRAC, 
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the  necessity  of  passing  all  these  functions  to  CALC  is  eliminated.  The  derivation 
of  these  direction  cosines  is  given  in  Appendix  VII. 

CONTROL  FUNCTIONS 

In  addition  to  the  orbital  calculations  just  described,  TRAC  subroutine 
functions  as  the  primary  controller  of  the  over-all  program.  The  control  constants 
for  the  program  are  contained  on  input  card  #2  which  is  read  into  the  program 
in  this  subroutine.  As  described  in  Section  II,  card  #2  contains  the  sixteen 
control  constants,  NUM  (1)  through  NUM  (16).  By  proper  choice  of  these 
constants  it  is  possible  to  control  the  integration  interval,  printout  intervals 
and  the  printout  increments  within  the  intervals,  and  the  selection  of  the  sub¬ 
routine  CALC  and  program  time  stop  options.  Assigned  control  constants  not 
used  in  TRAC  are  passed  on  to  the  appropriate  subroutines  called  by  TRAC. 

Some  of  the  storage  for  the  control  constants  is  not  utilized;  it  was  provided  to 
accommodate  future  program  modifications. 

NUM  (1)  is  the  number  of  ground  stations  in  the  program  and  controls  the 
loops  which  involve  station  parameter  calculations.  The  program  can  accom¬ 
modate  up  to  150  stations. 

NUM  (2)  equals  the  number  of  integration  increments  in  a  printout  increment. 
The  program  is  set  up  so  that  the  data  need  not  be  printed  out  after  every 
integration  interval.  For  example,  if  the  integration  increment,  which  is 
NUM  (8),  is  set  at  60  (in  seconds),  but  data  printout  is  desired  only  every  600 
seconds,  NUM  (2)  will  be  set  equal  to  600/60  =  10.  If  data  is  desired  at  every 
integration  increment,  NUM  (2)  will  be  set  equal  to  60/60  =  1. 

NUM  (3)  is  the  length  of  the  initial  printout  interval  for  the  variable  print¬ 
out  interval  option.  For  example,  it  may  be  desirable  to  print  out  data  every 
10  seconds  for  the  first  100  seconds  of  the  trajectory.  In  that  case,  NLTVI  (3) 
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will  be  set  equal  to  100.  The  variable  printout  interval  option  was  incorporated 
into  the  program  to  accommodate  the  high-loft,  suborbital  type  of  trajectory. 

In  this  case,  the  initial  and  final  parts  of  the  trajectory  are  changing  rapidly 
and  a  fine  printout  interval  is  desirable;  whereas,  in  the  middle  part  of  the  tra¬ 
jectory,  a  course  printout  is  usually  acceptable.  The  length  of  the  initial  phase 
is  controlled  by  NUM  (3).  The  altitude  at  the  time  in  trajectory  controlled  by 
NUM  (3)  is  then  stored.  When  the  trajectory  reaches  the  same  altitude  on  the 
return  path,  the  program  automatically  reverts  to  NUM  (3)  control.  Data  print¬ 
out  then  occurs  at  the  same  increment  established  for  the  initial  phase,  i.e.  , 

NUM  (2)  fixes  the  printout  increment  within  the  time  interval  established  by 
NUM  (3).  The  printout  option  will  be  bypassed  when  only  orbital  trajectories 
are  considered.  In  that  case,  the  trajectory  is  not  divided  into  intervals.  Print¬ 
out  occurs  at  the  printout  increment  established  by  NUM  (2)  until  the  program 
time  stop  is  reached. 

NUM  (4)  is  the  number  of  integration  increments  in  the  intermediate 
printout  interval.  It  does  not  affect  the  program  if  the  variable  printout  interval 
option  (controlled  by  NUM  (6) )  is  not  used. 

NUM  (5)  controls  the  final  program  time.  It  only  functions  as  a  program 
control  if  NUM  (6)  is  greater  than  one,  i.e.  ,  when  the  variable  printout  option 
is  by-passed.  If  NUM  (6)  is  zero  or  negative,  NUM  (5)  does  not  influence 
program  control.  The  product  NUM  (5)  x  NUM  (16),  executed  internally,  controls 
the  final  program  running  time.  For  example,  if  NUM  (16)  =  1,  NUM  (5) 
represents  the  number  of  seconds  to  final  orbit  time.  If  NUM  (16)  =,  60,  NUM  (5) 
represents  the  number  of  minutes  to  final  orbit  time. 

NUM  (6)  controls  the  selection  of  the  variable  printout  interval  option.  If 
the  variable  printout  option  is  by-passed,  the  program  will  run  until  time  equals 
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NUM  (5)  X  NUM  (16),  and  printout  will  occur  at  a  constant  increment  throughout 
the  program. 

NUM  (7)  controls  the  calling  of  the  subroutine  CALC  option.  If  it  is  greater 
than  zero,  the  CALC  option  is  selected. 

NUM  (8)  is  the  integration  interval  in  seconds. 

NUM  (9)  controls  the  calling  of  the  subroutine  DEBE  option.  If  it  is  greater 
than  zero,  the  DEBE  option  is  selected. 

NUM  (10)  through  NUM  (15)  are  unassigned. 

NUM  (16)  controls  the  final  program  time  as  explained  under  NUM  (5). 
Essentially,  NUM  (16)  permits  longer  program  running  times  without  upsetting 
the  data  card  input  format  established  for  the  NUM  (I)  series. 

Subroutine  TRAC  does  not  output  any  data  to  the  calling  program.  However, 
it  does  output  orbit  data  if  the  CALC  option  is  not  selected.  The  input  data 
required  from  the  calling  program  are  the  burnout  conditions  as  described  in 
Section  II.  These  data  are  normally  supplied  either  by  program  MAIN  or  the 
powered  flight  trajectory  simulation. 
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SECTION  IV 


GEOMETRICAL  COVERAGE  SIMULATION  -  SUBROUTINE  CALC 


The  primary  functions  of  subroutine  CALC  are  the  computation  of  the 
geometrical  relations  between  the  instantaneous  position  of  the  orbiting  vt  hide 
relative  to  the  ground  tracking  stations,  and  the  branch  control  of  the  DEBE 
subroutine  option  which  calculates  the  space  losses  for  various  vehicle-ground 
antenna  orientations.  Figure  8  is  a  flow  chart  illustrating  these  functions. 

STATION-VEHICLE  GEOMETRY 

In  subroutine  TRAC,  the  time-varying  vehicle  position  relative  to  the 
equatorial  plane  and  a  reference  meridian  (Greenwich)  was  established.  Latitude, 
longitude,  and  altitude  above  the  earth  were  the  three  basic  quantities  which 
completely  define  this  position.  In  subroutine  CALC  this  basic  position  data  is 
utilized  to  obtain  the  pertinent  quantities  relating  the  vehicle  position  to  a 
ground  tracking  station.  These  quantities  are  the  range  (RHO)  and  range  angle 
(PSI)  from  station  to  vehicle  and  the  azimuth  (AZD)  and  elevation  (ELE)  of  the 
line  of  sight  from  station  to  vehicle. 

From  this  data  and  a  knowledge  of  the  physical  limitation  of  the  tracking 
station  equipment,  tracking  coverage  of  a  vehicle  by  a  given  station  for  a  given 
ground  track  can  be  obtained  as  a  function  of  time.  Stations  below  the  minimum 
elevation  angle  are  rejected  by  the  program.  If  there  are  no  local  terrain 
limitations  or  inherent  physical  limitations  in  the  equipment,  the  minimum 
elevation  angle  (ELEM)  is  defined  as  the  angle  between  the  local  tangent  plane 
and  the  line  of  sight  tangent  to  the  earth' s  surface.  For  a  station  located  at 
sea  level,  the  minimum  elevation  is  zero.  However,  a  minimum  angle  of 


47 


ENTER  SUBROUTINE 

± 


subroutine 

STATEMENT  CARD 
(TRANSFER  ORBITAL 
POSITION  DATA  FROM  TRAC) 


DIMENSION  STATEMENT  CARD 
FOR 

ARRAY  STORAGE  ALLOCATION 


CODE 

^  PRINT  ^ 


INSTRUCTIONS 


INPUT  DATA 
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CALLING 
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Fig.  8  Flow  Diagram-Subroutine  CALC 
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5  degrees  is  often  used  to  allow  for  local  obstructions  in  the  case  of  the  optical 
line  of  sight  and  atmospheric  effects  in  the  case  of  the  electromagnetic  line  of 
sight.  In  any  event,  this  angle  is  selected  by  the  programmer  to  suit  the  appli¬ 
cation  and  is  inserted  into  the  program  along  with  other  station  data  (see  Table  5). 
The  nomenclature  of  the  parameters  defining  the  position  of  the  vehicle  relative 
to  the  earth  and  an  arbitrary  tracking  station  is  summarized  in  Fig.  9. 

Line-of-Sight  Range  from  Station  to  Vehicle  (RHO) 

The  cosine  law  applied  to  the  plane  triangle  OBS  in  Fig.  9  yields  the  range 
directly. 


RHO  =  f  (ALT)^  +  (RE  +  ALTS)^  -  2(ALT)  (RE  +  ALTS)  cos  (PSI)  (24) 

Range  Angle  (PSI) 

The  range  angle  is  defined  as  the  angle  between  the  radius  vectors  from 
the  center  of  the  earth  to  the  station  and  vehicle,  respectively.  It  can  vary 
between  0  and  180  degrees.  Application  of  the  cosine  law  to  the  spherical 
triangle  PSB  in  Fig.  9  yields  (actually,  only  the  cosine  of  PSI  and  not  PSI  itself 
is  required  for  the  simulation) 

cos  (PSI)  =  cos  (90  -  ALAT)  cos  (90  -  RLAT)  +  sin  (90  -  ALAT) 

X  sin  (90  -  RLAT)  cos  (RLONG  -  ALONG)  (25) 

Substituting  compliments  and  expanding  the  cos  (RLONG  -  ALONG)  gives; 

(JOS  (PSI)  =  sin  (ALAT)  sin  (RLAT)  +  cos  (ALAT)  cos  (RLAT)  [cos  (ALONG) 
cos  (RLONG)  -  sin  (ALONG)  sin  (RLONG)]  (26) 
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ALONG 


STATION 


RLAT 


-RUONG 


ASNOa 


ORBIT 


VEHICLE 


Si  =  STATION  LOCATION 

B  =  VEHICLE  LOCATION 

ALONG:  LONGITUDE  OF  STATION 
RLONG=  LONGITUDE  OF  VEHICLE 
ALAT  =  LATITUDE  OF  STATION 
RLAT  =  LATITUDE  OF  VEHICLE 

AZD  =  AZIMUTH  OF  VEHICLE  FROM  STATION  MEASURED  FROM  NORTH  CLOCKWISE 
PSI  =  RANGE  ANGLE  OR  ANGLE  BETWEEN  STATION  AND  VEHICLE  MEASURED  FROM 
EARTH  CENTER 

RHO  =  DISTANCE  BETWEEN  STATION  AND  VEHICLE 
ALT  =  ALTITUDE  FROM  EARTH  CENTER  TO  VEHICLE 
ALTS  :  ALTITUDE  FROM  EARTH  SURFACE  TO  SENSOR 

RE  =  RADIUS  OF  THE  EARTH 

OSi  =  RE  +  ALTS 


Fig.  9  Station-Vehicle  Geometry 
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In  Fortran  notation  this  becomes 


CPSI  =  SLA(L)  *  SLAR  +  CLA(L)  *  CLAR  *  (CLO(L)  *  CLOR  +  SLO(L)  *  SLOR) 

(27) 

where  L  is  the  station  index. 

Azimuth  Angle  (AZD) 

The  application  of  the  cosine  law  to  the  spherical  triangle  PSB  in  Fig.  9 

yields 


cos  (AZR) 


sin  (RLAT)  -  cos  (PSI)  sin  (ALAT) 
sin  (PSI)  cos  (ALAT) 


(28) 


In  Fortran  cos  (AZR)  =  PHI,  and  f  1  -  cos^  (PSI)  is  substituted  for  sin  (PSI) 
giving 


PHI  -  (SLAR  -  CPSI  *  SLA  (L))/CLA(L)  *  SQRT  (1.0  -CPSI**2)  (29) 

.'.  AZR  =  A  cos  (PHI) 

and  AZD  -  AZR  *  CONVI  (30) 

where  CONVI  is  a  conversion  factor  from  radius  to  degrees. 

Because  the  azimuth  or  heading  angle  between  the  tracking  station  and  the 
vehicle  in  orbit  is  measured  from  north  in  a  clockwise  direction,  a  discontinuity 
in  the  heading  angle  can  occur  for  certain  orbital  paths  relative  to  the  ground 
station.  If  the  vehicle  latitude,  for  example,  is  greater  than  the  station  latitude, 
and  the  orbital  path  remains  north  of  the  station,  the  measured  azimuth  track 
angle  will  vary  from  some  value  between  +  90  and  0  degrees  through  0  to  some 
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value  between  360  and  270  degrees.  The  appropriate  branch  logic  control  must 
be  provided  in  the  simulation  to  account  for  such  orbital  tracking  situations. 

The  permutations  in  the  heading  angle  which  may  occur  for  possible  '  hide 
tracking  station  locations  are  illustrated  in  Fig.  10.  Two  types  of  orbital  paths 
have  to  be  given  special  consideration.  In  the  first  type,  the  orbit  passes  from 
the  first  to  fourth  quadrant  in  the  heading  angle  but  does  not  cross  the  discon¬ 
tinuity  in  longitude  caused  by  measuring  longitude  east  and  west  of  the  Greenwich 
meridian.  In  the  second  type  the  orbit  passes  from  the  first  to  the  fourth  quadrant 
in  heading  and  also  crosses  the  discontinuity  in  the  longitude  at  ±  180  degrees. 

As  shown  in  Fig.  10,  the  proper  value  of  the  heading  angle  AZD  can  be 
obtained  by  branching  on  tests  which  compare  the  longitude  of  the  station  and 
vehicle.  The  following  branching  logic  takes  into  account  all  possible  permutations. 
Let  be  the  azimuth  angle  computed  by  Eq.  28,  then  if 


RLONG  -  ALONG 


RLONG  -  ALONG 


RLONG 


RLONG 

RLONG 


RLONG 


(ALONG  +  180) 

(ALONG  +  180) 

(ALONG  -  180) 

(ALONG  -  180) 


0  —  AZD  =  /3 

0  -*  AZD=  360  -  ^ 

0  -*AZD  =13 

0  —AZD  =  360  -  ^ 


The  Fortran  listing  for  this  portion  of  the  program  is  given  on  lines  Cl 6  through 
C23  in  Appendix  HI. 
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a.  VEHICLE  LONGITUDE  EAST  OF  GREENWICH 
10  <  RLONG  <  180) 


0<  ALONG  <RLONG 


IF 

RLONG<ALONG<  160 


OR 


OR 


RLONG-(ALONG  +180X0 
RLONG  -  ALONG  >0 


THEN 


-  (180-  RLONGX  ALONG  <0 


OR 


RLONG- (ALONG  +  180X0 
RLONG- ALONG  >0 

THEN 


AZD=  /3 


IF 

-I80<AL0N6<-(I80-RL0NG) 

OR 

' 

RLONG -(ALONG +180)  >  0 
RLONG  -  ALONG  >0 

THEN 

AZD=360-)9 


Fig.  10  Resolution  of  Heading  Ambiguity 
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b.  VEHICLE  LONGITUDE  WEST  OF  GREENWICH 
-I80<RL0NG  <  0 


0<AL0NG<  (I80  +  RL0NG) 


OR 

RLONG  -  (ALONG - 1 80) >  0 
RLONG  -  ALONG  <  0 


RLONG<ALONG<  0 
OR 


RLONG-(ALONG  +  180X0 
RLONG- ALONG  <0 


THEN 

AZD*360-i8 


180+ RLONG  <  ALONG  <  180 
OR 

RLONG -(ALONG -180) <  0 
RLONG  -  ALONG  <  0 


THEN  -O 


-  I80<AL0NG<  RLONG 
OR 


RLONG  -( ALONG  +180)  >  0 
RLONG -ALONG  >0 


THEN 
AZD  =  i8 


Fig.  10  (Continued)  Resolution  of  Heading  Ambiguity 
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Elevation  Angle  (ELE) 


The  elevation  angle  can  be  found  by  applying  the  law  of  sines  to  the  plane 
triangle  SBS'  in  Fig.  11  (which  is  an  expanded  detail  of  the  triangle  OSB  of 
Fig.  9). 


sin  (ELE) 
BS' 


sin  p 

BS' 


sin  p 
RHO 


sin  (180  -  (90  -  PSI))  =  cos  PSI 


(RE  +  ALTS) 
cos  PSI 


(31) 

(32) 

(33) 
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Substituting  Eq.  (32)  and  (33)  into  (31)  and  solving  for  sin  (ELE)  gives 


sin  (ELE) 


ALT  cos  (PSD  -  (RE  +  ALTS) 

RHO 


(34) 


In  Fortran  notation 


SELE  =  [ALT  *  CPSI  -  (RE  +  ALTS)  ]  /RHO 


(35) 


ELE  =  ASIN  (SELE) 


(36) 


ANTENNA  ASPECT  ANGLES  (ANGLE,  FNGLE,  ALPHA) 

In  addition  to  vehicle- station  geometry,  CALC  subroutine  calculates  the 
following  angles  which  define  various  vehicle-sensor  antenna  axes  orientations. 

ANGLE  -  the  angle  between  the  line  of  sight  from  station  to  vehicle  and 
the  axis  of  the  spin- stabilized  vehicle  antenna. 

FNGLE  -  the  angle  between  the  line  of  sight  from  station  to  vehicle  and 
the  axis  of  an  earth-center-pointing  vehicle  antenna. 

ALPHA  -  the  angle  between  the  line  of  sight  from  station  to  vehicle  and 
the  fixed  axis  of  the  i-th  station. 

The  cosine  of  each  of  these  angles  is  found  by  taking  the  dot  product  of 
unit  vectors  along  the  line  of  sight  and  the  respective  antenna  axes.  The  required 
direction  cosines  which  define  these  unit  vectors  are  derived  in  Appendix  VII. 

The  above  angles  are  given  by  the  inverse  cosine  of  Eqs.  (73),  (74),  and  (75)  as 
follows. 
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ANGLE  =  ACOS  (XI  *  XRHO  +  YI  *  YRHO  +  Zl  *  ZRHO) 


(37) 


FNGLE  =  ACOS  (XF  *  XRHO  +  YF  *  YRHO  +  ZF  *  ZRHO)  (38) 

ALPHA  =  ACOS  (XAN(L)  *  XRHO  +  YAN  (L)  *  YRHO  +  ZAN  (L)  *  ZRHO)  (39) 

where 


XI,  YI,  ZI 

XF,  YF,  ZF 


XAN(L),  YAN(L),  ZAN(L) 


RATES  OF  CHANGE 


=  the  direction  cosines  of  the  spin- stabilized 
antenna 

=  the  direction  cosines  of  the  spaceborne 

antenna  pointing  toward  the  center  of  the  earth 

=  the  direction  cosines  of  the  fixed  antenna 
axis  of  the  L-th  ground  station  (or  sensor) 


It  is  often  desirable  to  know  how  fast  a  variable  is  changing.  Therefore, 
the  time  rates  of  change  of  RHO,  AZD,  ELE,  ANGLE,  FNGLE,  and  ALPHA  are 
computed  and  printed  out.  In  this  program  a  simple  linear  difference  between 
computation  points  is  used.  It  is  planned  to  incorporate  expressions  for  the 
instantaneous  rates  of  change  of  these  quantities  into  the  program  in  the  very 
near  future. 
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SECTION  V 


ELECTROMAGNETIC  COVERAGE  SIMULATION  -  SUBROUTINE  DEBE 


All  antenna  gain  and  space-loss  calculations  are  made  in  one  subroutine 
called  DEBE.  The  use  of  this  subroutine  is  optional  and  is  controlled  by  the 
fixed  point  number,  NUM(9),  on  the  control  parameter  input  card.  In  its  pres¬ 
ent  form,  the  subroutine  calculates  the  total  losses  for  eight  different  vehicle- 
ground  station  antenna  configurations.  However,  by  consolidating  all  of  the 
antenna  calculations  in  one  subroutine,  program  flexibility  is  achieved,  and 
other  antenna  orientations  and  performances  can  be  studied  by  simply  inserting 
the  appropriate  DEBE  subroutine.  A  fimctional  flow  diagram  of  the  required  t 
calculations  is  shown  in  Fig.  12. 

The  total  losses  are  made  up  of  the  space  loss  (SLOS)  and  the  gains  (losses) 
of  the  various  antennas.  Free-space  conditions  are  assumed;  that  is,  interfer¬ 
ence  and  diffraction  phenomenon  due  to  meterologlcal  conditions  and  other  fac¬ 
tors  are  neglected.  The  various  antenna  combinations  are  summarized  in 
Table  11,  and  the  associated  geometry  is  shown  in  Fig.  13. 


SPACE  LOSS  (SLOS) 


The  total  signal  p>ower  received  is  given  by  the  standard  one-way,  free- 
space  propagation  equation  for  the  ratio  of  p>ower  received  to  power  trans¬ 
mitted.^^’ 


P 

r 


grgt 


(40) 
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ENTER  SUBROUTINE 
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Functional  Flow  Diagram— Subroutine  DEBE  Fig.  13  Geometry  for  Space-Loss  Calculati 


where 


P^,  =  received  power  and  transmitted  power,  respectively. 

X  =  wavelength,  in  the  same  units  as  R. 

g^,  g^  =  power  gains  of  the  receiving  and  transmitting  antennas, 
respectively. 

R  =  line  of  sight  range. 

Table  11 

Antenna  Combinations 


Associated 

Vehicle  Antenna 

Ground  Antenna 

Total 

Loss* 

Direction 
(See  Fig.  13) 

Gain 

Direction 
(See  Fig.  13) 

Gain 

DBl 

S.A. 

GANGLE 

F.A. 

GALPHA 

DB2 

S.  A. 

GANGLE 

RHO 

DBA  (L) 

DB3 

E.  C. 

GFNGLE 

F.A. 

GALPHA 

DB4 

E.  C. 

GFNGLE 

'  RHO 

DBA  (L) 

DBS 

RHO 

DBG 

F.A. 

GALPHA 

DB6 

RHO 

DBG 

RHO 

DBA  (L) 

DB7 

S.A. 

GANGLE 

Isotropic 

0 

DBS 

E.  C. 

GFNGLE 

Isotropic 

0 

*Total  loss  equal  space  loss  plus  the  gains  of  the  antenna  pair,  i.  e. , 


DBl  =  SLOS  +  GANGLE  +  GALPHA. 


61 


For  purposes  of  computational  flexibility  the  total  loss  as  given  by  Eq.  (40) 
is  computed  term-by-term  in  decibels  (10  log  P^/P^)  and  added  to  find  the 
total  loss.  The  total  loss  for  unity  antenna  gains  is  designated  as  the  spa-,  loss 
(SLOS)  in  the  simulation,  thus 


SLOS  =  (j\\ 

\47r  R/ 

where  the  wave  length  (X )  is  related  to  the  system  frequency  (FREQ)  by  the 

speed  of  light  (c)  in  vacuo.  If  the  range  (RHO)  is  given  in  nautical  miles  and 

FREQ  in  megacycles,  then  the  speed  of  light  which  is  299,  776  x  10^  meters/sec. 

— 6 

must  be  expressed  in  nautical  miles  per  second  times  10  or 


=  /— 

^FREQ 


FREQ  Att  RHO 


)■ 


(41) 


c 


(299.  776)  (3.28083)  /n.  m. 
6076.115  \sec. 


X 


) 


In  the  simulation,  space  loss  is  computed  as  the  sum  of  two  components, 
the  frequency-dependent  space  loss  (FDSL),  and  the  range-dependent  space 
loss  (RDSL)  expressed  in  decibels.  Thus, 


SLOS  =  FDSL  +  RDSL 


where 


FDSL  =  20.  *ALOG10[  (299.776  *  3.  28083/6076. 115)/FREQ]  (42) 
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and 


RDSL  =  -20.  *  ALOGIO  (4.  *  3.  1415927  *  RHO)  (43) 

ANTENNA  GAINS 
On-axis  Gains  (DBG,  DBA) 

For  antennas  which  radiate  (or  receive)  a  well-defined  beam,  the  funda¬ 
mental  relation  which  connects  the  on-axis  or  maximum  gain  G^,  the  area  of 
the  aperture  A,  and  the  wavelength  is: 


G  =  f 
o 


47r  A 


(44) 


The  dimensionless  factor  f  is  equal  to  1  if  the  excitation  is  uniform  to  phase 
and  intensity  over  the  whole  aperture;  in  actual  antennas,  f  is  oten  as  large 
as  0.  8  and  is  rarely  less  than  0.  5.  The  connection  between  gain  and  beamwidth 
is  as  follows.  Using  an  aperture  of  dimensions  d  in  both  directions,  a  beam 
may  be  formed  whose  angular  width*,  determined  by  diffraction,  is  very  nearly 

BW  =  -^ 
d 

2 

If  the  aperture  is  square  A  =  d  ,  thus 
the  gain  becomes 

G  =  f 
o 


radians 


(45} 


substituting  for  X  and  A  in  Eq.  (44), 


47r 


(BW) 


(46) 


*The  beamwidth  is  the  angular  interval  between  two  directions  for  which  the  gain 
is  half  of  the  maximum  gain. 
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TT  2 

If  the  aperture  is  circular  A  =  ~  d  ,  the  gain  would  be 


2 


G  =  f 


(47) 


o 


In  DEBE  subroutine,  Eq.  (46)  is  used  for  the  maximum  gain  of  the  transmitting 
and  receiving  antennas.  The  transmitting  antenna  is  considered  to  be  on  the 
space  vehicle,  and  since  it  is  a  transmitting  antenna  close  to  the  source  of 
illumination,  a  value  of  0.  8  was  selected  for  the  factor  f.  This  reasoning  is 
also  consistent  with  the  type  of  antenna  likely  to  be  used  on  the  vehicle,  such 
as  a  helix  or  a  dipole,  which  have  higher  efficiencies  than  a  dish-type  likely  to 
be  used  on  the  ground.  Thus,  the  maximum  gain  of  the  vehicle  antenna  is 


/  3.5449l\^ 

\  BW  / 


(48) 


For  the  receiving  antenna  on  the  ground  where  the  illumination  is  likely  to  be 
a  less  uniform  factor,  f  =  0.  6  was  selected.  The  ground- antenna  maximum 
gain  is,  therefore. 


(49) 


where  in  both  instances  BW  is  the  beamwidth  between  half-power  points  in 
radians.  Of  course,  for  the  isotropic-ground-antenna  case,  the  gain  is  unity 
in  all  directions  or  zero  db. 
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Off-axis  Gains  (GANGLE,  GFNGLE,  GALPHA) 


Cases  are  considered  in  which  the  vehicle  and  ground  antennas  are  not 
completely  directional  relative  to  each  other.  To  handle  these  cases,  it  is 
necessary  to  simulate  the  antenna  pattern.  There  are  an  unlimited  number  of 
possible  patterns.  However,  one  that  is  fairly  representative  of  many  antennas 
and  also  was  appropriate  to  the  original  problem  (see  Section  2,  General 
Description)  is  the  uniform  distribution,  which,  in  terms  of  the  power  gain,  is 
given  by 


where  a  is  the  angle  off  the  axis.  It  can  be  any  of  the  angles  (ANGLE,  FNGLE, 
ALPHA)  shown  in  Fig.  13.  The  half-power  beamwidth  can  be  expressed  as 


BW 


K 


\ 

d 


which,  when  substituted  in  Eq.  (50),  gives 


(51) 


(52) 
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The  constant,  K,  can  be  determined  from  the  condition  at  the  half-power  point, 


i.  e. ,  when  o; 


1 

2 


BW, 


since  the  factor  multiplying 


G 

o 


must  be  equal 


Equation  (53)  is  satisfied  by  a  value 


Ktt  =  2.78332 


Thus,  the  final  expression  becomes 


G 


G 


sm 


/2.  78332a  \ 

V  BW  / 


2.  783320; 


BW 


In  Fortran  notation,  the  three  antenna  gains  are 
the  spin-stabilized  antenna  gain 

GANGLE  =  10.  ♦  ALOGIO  (CON2  *  ARG2) 
the  fixed-ground  antenna  gain 

GALPHA  =  10.  *  ALOGIO  (CONI  *  ARG) 
the  earth- centered  antenna  gain 

GFNGLE  =  10.  *  ALOGIO  (CON2  *  ARG3) 


1 

to  -  . 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 


66 


where 


CONI  =  DBA  (L)  =  (3.54491/BW  (L)  *  CONV))  **  2  *  0.  6 

CON2  =  DBG  =  (3.54491/BWl)  **2*0.8 

ARG  =  (SIN  (AN1)/AN1)  **  2 
ARG2  =  (SIN  (AN2)/AN2)  **  2 
ARG3  =  (SIN  (AN3/AN3)  **  2 

where 

ANl  =  2.78332  *  ALPHA/BW  (L) 

AN2  =  2.  78332  *  ANGLE/BWG 

AN3  =  2.78332  *  FNGLE/BWG 

CONV  is  a  conversion  factor  from  degrees  to  radians 
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SECTION  VI 


SOLUTION  OF  THE  ORBITAL  EQUATION  -  SUBROUTINE  RUNG 


Subroutine  RUNG  performs  the  numerical  integration  of  the  Keplerian 
equation  for  the  rate  of  change  of  the  true  anomaly  of  an  elliptical  orbit.  The 
equation  in  Fortran  notation  is 

[1.  +  ECC  X  cos  (VEE)  ]  ^  =  F  (VEE)  ( 

P 


where 


VEE  =  true  anomaly 
ECC  =  eccentricity 

3  -2 

GMU  =  gravitational  mass  unit  (nm  sec  ) 

P  =  semi-latus  rectum 

For  a  given  elliptical  orbit  about  a  spherical  earth,  the  vehicle  position 
in  the  orbit  is  completely  defined  by  the  true  anomaly  (VEE)  which  is  the  angular 
position  measured  from  the  line  of  perigee.  Hence,  the  numerical  integration 
procedure  which  calculates  the  true  anomaly  establishes  the  position  of  the 
vehicle  in  a  given  elliptical  orbit  at  the  next  time  increment. 

The  input  data  to  the  subroutine  is  the  true  anomaly  at  the  present  instant 
of  tinio  (VEE(K)  and  the  orbital  parameters,  K,  P,  and  ECC.  The  subroutine 


69 


then  calculates  the  true  anomaly  VEE  (K+1)  at  the  next  instant  of  time  as 
determined  by  the  integration  interval  which  can  be  arbitrarily  selected. 

The  subroutine  integrates  the  true  anomaly  equation  given  above  ’  ’ 
following  fourth-order  Runge-Kutta  integration  technique. 

VEE  (K+1)  =  VEE(K)  +  (vEER(l)  +  2  (VEER  (2)  +  VEER  (3))  +  VEER(4))/6  (60) 

where 


VEE  (K) 

the  value  of  VEE  at  the  start  of  the  time  increment 

VEE  (K  +  1)  = 

the  value  of  VEE  at  the  end  of  the  time  increment 

VEER(l)  = 

F  (VEE  (K))  *  TNUM 

VEER  (2)  = 

F,VEE(K)  + 

*  TNUM 

VEER  (3)  = 

VEER  (2)  \ 
F  (VEE  (K)  +  2 

*  TNUM 

VEER  (4)  = 

F(VEE(K)  +  VEER  (3)) 

*  TNUM 

TNUM 

integration  interval 

The  calculated  value  of  VEE  (K  +  1)  is  passed  to  subroutine  TRAC  where 
the  calculations  involving  true  anomaly  are  made.  After  these  calculations, 

VEE  (K  +  1)  is  returned  to  this  subroutine  and  becomes  VEE  (K)  in  the  calculation 
of  the  true  anomaly  at  the  succeeding  time  interval.  This  cycling  is  repeated 
until  the  program  stops. 

The  errors  in  the  true  anomaly  calculated  by  this  numerical-integration 
technique  have  not  been  analytically  studied  in  detail.  However,  the  equation 
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for  the  rate  of  change  of  the  true  anomaly  is  well  behaved,  and  the  integration 
technique  is  adequate  for  the  engineering  applications  for  which  the  program  is 
intended.  A  functional  flow  diagram  of  the  subroutine  is  given  in  Fig.  14. 


L.  E.  Wilkie 


4 


J.  W.  Stevens 
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ENTER  SUBROUTINE 


RETURN  TO  CALLING 
PROGRAM 


.  14  Functional  Flow  Diagram — RUNGE-KUTTA  Subroutine 
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APPENDIX  I 


FORTRAN  IV  LISTING 
PROGRAM  MAIN 


T  SU0TYPE«FORTRAN«LMAP 

- C“  -  WOOTFrEO  SPACfc"  VEHICLE  TPACkInG  AND  SENSOR  COVERAGE  PROGRAM 

C  BY  J,STeVENS»L«WILKIE*D23*MITRE  CORPOR AT I  ON • BEDFORD • MASS 

CMAIN 

C  THIS  PROGRAM  IS  MADE  UP  OF  THE  FOLLOWING  SUBROuTInES- 
C  TRAC  -COMPUTES  THE  ORBITAL  POSITION  AND  ATT  I TuDE  OF  THE  VEHICLE. 

C  RUNG  -COMPUTES  THE  TRUE  ANOMALY  AS  A  FUNCTION  OF  T I  ME  _  _ _ _ 

"C"  “  By  a  fourth 'order  runge  kutta  numerical  Tntegr'aticn  method 

CALC  -COMPUTES  THE  GEOMETRICAL  RELATIONS  BETWEEN  THE  VEHICLE 
AND  THE  GROUND  SENSORS. 

DEBE  -COMPUTES  THE  SPACE  LOSSES  AND  THE  GAINS  OF  VARIOUS  VEHICLE- 


GROUND  STATION  ANTENNA  CONFIGURATIONS. 

1  READ  2.HBO. VBOtGBO.BOHEADrBOLAT .BOLONl tBIG.BlH  M  I. 

2  FORTOTlSFld.f)  ”  '  M2. 

CALL  TRAC  ( HBO ♦ VBO . GBO ♦ BOHE AD . BOLAT . BOLON I . B I G • B I H )  M  3. 

IF(VBO)3.1.1  M  4. 

3  RETURN  M  5* 

END  M  6. 
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APPENDIX  II 


FORTRAN  IV  LISTING 


SUBROUTINE  TRAC 


r 

C 

SYMBOL 

_ LIS 

UNITS 

IN  INF  )  HA^  1 _ 

OEF INI T ION 

C 

alt 

(N  M  I  ) 

altitude  FROM  THE  CENTER  OF  ThE  EARTH 

c 

alt  I 

(N  MI  ) 

altituoe  from  the  Surface  of  the  earth 

c 

ASNOOE 

(  RAD  ) 

LONGITUDE  of  the  ASENDING  NODE • ME ASURED  FROM  THE 

c 

GREENWICH  MERIDIAN  FOR  A  NON-ROTATING  EARTH 

c 

BETA 

(DEG) 

ARGUMENT  OF  .PERIOEE   

c 

BIG 

(OEG  ) 

BURNOUT  SPIN  AXIS  ATT  I  TUOE  •  POS  I  T  I  VE  UP  FROM  LOC;  L. 

c 

HORIZONTAL  PLANE 

c 

BIH 

(OEG) 

BURNOUT  SPIN  AXIS  HE AO  I NG t PO S I T I VE  CLOCKWISE  FROM 

c 

NORTH 

c 

BOHEAO 

(OEG  ) 

BURNOUT  HEADING. POSIT  I VE  CLOCKWISE  FROM  NORTH 

c 

BOLAT 

(OEG) 

BURNOUT  LATITUDE.POSITIVE  NORTH  EROM  THEELQUATOR 

c 

BOLONG 

(  OEG  ) 

burnout  long  I TUDE .POSI T I VE  EAST  FROM  GREENWICH 

c 

Cl 

( NQr^  ) 

CONSTANT  FACTOR  OF  THE  SP IN-STAB I  LI 2EX)  ANTENNA 

c 

DIRECTION  COSINES 

c 

C2 

( NONE ) 

CONSTANT  factor  OF  THE  SPIN-STABILIZED  ANTENNA 

c 

direction  cosines 

c 

ECC 

(NONE ) 

ECCENTRICITY  OF  THE  ORBIT 

c 

GBO 

(OEG) 

BURNOUT  Plight  path  angle .positive  up  from  local 

c 

HOR IZOnTALPLANE _  _ 

c 

HBO 

(N  M  I  ) 

BURNOUT  alt  I TUDE. positive  UP  FROM  SURFACE  OF  FAPTH 

c 

NUM (  1  ) 

( NONE ) 

NUMBER  OF  GROUND _SE NS OR S 

c 

NUM ( 2 ) 

NUMBER  OF  INTEGRATIONS  IN  THE  INITIAL  ANO 

c 

final  PRINTOUT  I NTERVAl 

c 

NUM ( 3 ) 

(SEC) 

ENO  OF  FIRST  PRINTOUT  INTERVAL 

c 

NUM ( 4  ) 

NUMBER  OF  integrations  IN  THE  INTERMEDIATE 

c 

PRINTOUT  interval 

c 

NUM ( 5  ) 

(  SEC  ) 

final  time 

c 

NUM ( 6  ) 

( NONE ) 

CONTROL  VARIABLE  TO  SELECT  TImE  STOP  OPTION 

c 

NUM ( 7 ) 

( NONE ) 

CONTROL  VARIABLE  TO  SELECT  CAl  C 1  OPTION 

c 

NUM ( 8 ) 

integration  interval  in  seconds 

c 

NUM ( 9 ) 

(NONE ) 

CONTROL  VARIABLE  To  SELECT  OEpE I  OPTION 

c 

NUM (  I  6 ) 

(NONE ) 

multiplication  factor  to  increase  NUM(5).IT 

c 

MUST  BE  A_P05ITIVF  INTEGER 

c 

NUM  2 

(SEC) 

INITIAL  PRINTOUT  INTERVAL 

c 

NUM4 

(SEC  ) 

INTERMEOIATE  PRINTOUT  INTERVAL 

c 

omega 

(RAO  ) 

argument  of  perigee 

c 

P 

(N  M  I  ) 

SEMI-LATUS  RECTUM 

c 

PLA  IN 

(OEG) 

orbital  inclination 

c 

SNOOE 

(DEG) 

ASNOOE  IN  OEGREES 

c 

VBO 

(FT/SEC ) 

BURNOUT  VELOCITY 

c 

VEE 

(RAO  ) 

TRUE  ANOMALY 

c 

VEE  I 

(OEG  ) 

TRUE  ANOMALY 

c 

XI 

direction  COSINES  OF  THE  SP I N-ST AB I L I  ZED  ANJENNA 

c 

YI 

(NONE ) 

RELATIVE  TO  AN  X.Y.2.FRAME  CENTERED  IN  THE  EARTH. 

c 

21 

2  ANO  Y  ARE  POSITIVE  NORTH  AND  E A ST f RESPECT LVELY 

c 

this  Subroutine 

CALCULATES  VEHICLE  POSITION  IN  GEOCENTRIC  COORO I N-  T 

c 

ates.the 

I NPUTS 

ARE  The  burnout  conditions  WHICH  may  be  obtained  t 

C  ElTHrp  rpoM  MAIN 
C  1  he  SUBROO  r ! NC  U 
C  CALLING  STATEMEN 

subroutine  TRAC 


BOOST  PR0GRAM5( 20R3nFGRFFS 
TO  MAKE  IT  compatible  WITH 


OF  FRFFHOM  ) 

general 


OR  THE 
SES  TRAC 
T 

(HBO* VB0*GB0*B0HEAD*B0LAT •bOLON I •BIG»BIH  ) 
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DIMENSION  NUM{ 16 J . VFF 

T 

2  • 

POLONG=POLONI 

IF (  I  ST ART + I  )  I  91  •193^191 

T 

2.01 

I  91 

C 

ISTART=-I 

physical  constants 

^.21 

PIE  =  3.  I  4  I  5927 

T 

2  •  3 

HAPIF=  1  .5707963 

T 

2*4 

CONV=O.OI 745329 

T 

3. 

CONVI  =57.2957B _ _ 

T 

4  . 

CFEET  =  6076. I  I  5 

T 

5. 

CNM I =0*0001 6457BB 

T 

6  • 

RE  =  34  3B. I  49 

T 

7. 

WE=2.0*PI E/B6400. 

T 

A* 

GMU=62750.2I 

T 

9. 

.REA0_192*  (NUM(H  .U  =  l  f  J6) _  _ 

T 

9*5 

I  92 

format (1615) 

T 

9.6 

NUM2=NUM { 2 ) *NUM (B ) 

T 

IQ*. 

NUM4  =NUM { 4 ) *NUM { 8 ) 

T 

I  I  • 

NUM  (  5  1  =NUM  I  5  ^  *NUMf  .1  6  i 

T 

*  ?  T  O  1 

PRINT  69 

T 

I  I  •  I 

PRINT  B7  f  (  NUM  (L  )  •L=  I  f  )  _  _ 

T 

I  I  •  12 

I  93 

PR  I  NT  72  f  HBO • VBO • GBO • BOHEAO • BOLAT • BOLONG • b I G , B I H 

T 

12.  I 

HBO=HbO+RE 

T 

12.2 

C 

VBO=VBO*CNM I 

CONVERT  degrees  TO  RAO  I ANS 

T 

13. 

GBO=CONV*GBO 

T 

I  4 

BOHEAR=CONV*BOHEAO 

T 

15. 

B0LAT  =  C0NV*B0LAT  . .  ' 

T 

16. 

BOLONG=CONV*BOLONG 

T 

17* 

BIG=CONV*BIG 

T 

is 

BIH=C0NV*BIH 

T 

19 

C 

compute  Sin  and  cos  of  burnout  parameters 

CHEBO=COS ( BOHEAR ) 

T 

20  • 

SHEBO  =  S I N ( BOHE  A  R ) 

T 

21 

C0BIG=C0S(BIG> 

T 

2^  • 

SOB  I G  =  S I N ( B I G  > 

f 

23. 

COB  I H  =  COS(B I H ) 

_ I 

24. 

SOBIH=SIN(BIH) 

T 

25* 

COGBO=COS(GBO) 

T 

26. 

SOGBO=S IN (GBO ) 

T 

21. 

CLABO  =  COS ( BOLAT  ) 

T 

28. 

SLABO=SI N(BOLAT ) 

T 

29. 

CLOBO=COS(POLONG ) 

T 

30^ 

SLOBO  =  S I N ( POLONG ) 

T 

31  . 

C I =  SOB  I G*CLABO-COB I G*COB I H^SLABO 

a2* 

C2=C0BIG*S0BIH 

T 

33. 

ZI =SOBI G*SLABO+COB I G*COB I H^CLABO 

T 

33.1 

~c 

COMPUTE  ORBITAL  PARAMETER’S  «  ECCENT9  I  C  I  T  Y  •  I  NCL  I  NAT  I  ON  .  TRUE  ANOMALY 

A=HB0*VB0**2/GMU 

T 

34  . 

ECC=SORT ( ( (A- I • 0) *COGBO ) **2+ SOGBO** 2 ) 

T 

35. 

►  _  _ 

PLA [N=ACOS(CLABO*SHEBO ) *CONV I  _ 

T 

36. 

t 

D=A*C0GB0**2-I *0 

T 

37. 

\  _ 

VEE(  I  )=ATAN(A*SOGBPj»COGBP/_D) 

T 

38* 

IF(D)20I *202,202 

T 

39. 

201 

VEE( I )=Vre( I )+PIE 

T 

40. 

202 

C 

C 

C 

CONT I NUE 

COMPUTE  OMEGA  *  ARGUMENT  OF  PERIGEE  RELATIVE  TO  ASCENDING  NOOE 
COMPUTE  ASNODE* LONG  I TUDE  OF  ASCENO I NG  NODE  MEASURED  FROM  THE 
GREENWICH  MERIDIAN  ANO  RELATIVE  TO  A  NON-ROTATING  EARTH 

T 

4  I  • 

C I NL  =  CLABO*SHEBO 

T 

42  . 

76 


c 

80  1 
810 

SI NL=SORT ( 1 .-C I NL**2 ) 

slain=slapo/sinl 

Ab=ASIN (SLAIN ) 

AC=AS1N( SLAIN*5H580) 

continue 

USE  boheao  to  resolve  quadrant  ambiguity  of  Omega  and  asnode 

T 

T 

T 

T 

T 

43. 

44  . 

45  . 

52. 

53. 

1F(B0HEAD-90.0)8I I .81 1 .805 

T 

54. 

805 

IF (BONE AD- 1 80.0 ) 806. 806. 807 

T 

55. 

807 

IF (bOHEAD-27b.O ) 806 .806.81 1 

T 

56. 

81  1 

OMEGA=-VEE( 1 )+AB 

T 

57. 

ASNODE =BOLONG- AC 

T 

58. 

G0TO804 

T 

59. 

806 

OMEGA=-VEE( 1 )-AB+PlE 

T 

60. 

ASNODE=POLONG+AC+P IE 

T 

6  1. 

AB=P IF-AP 

T 

6  1.1 

804 

COnT I NUE 

T 

62. 

C0M=C0S(0MEGA) 

T 

63. 

SOM=S IN (OMEGA) 

T 

64. 

c 

COMPUTE  constant  FACTORS  IN  ORBITAL  POSITION  EOUATIONS 

C3=S0M*SINL 

T 

65. 

C4=C0M*SI NL 

T 

66  . 

P=HB0* ( 1 .+ECC*COS( VEE ( 1 ) ) ) 

T 

67, 

VEE I =VEE ( 1 ) *CONV I 

T 

68. 

B=SORT (GMU ) /P** 1 .5 

T 

69. 

BETA=C0NV I *OMEGA 

T 

70. 

SNODE  =CONV I  * ASNODE 

PR  I  nt^aa^eccTplaTnTp.  v^  1  .reta.snode” 

T 

71.1 

c 

COMPUTE  TIME  varying  FACTORS  IN  ORBITAL  POSITION  EOUATIONS. 

c 

LAT ITUDE. longitude. AND  ALTITUDE  OF  SPACE  VEHICLE  RELATIVE  TO 

c 

MOVING  earth 

1  =  1 

T 

71.1 

1TIME=0 

T 

72. 

700 

CVEE=C0S( VEE( 1 ) ) 

T 

73. 

SvEE=S1N( VEE ( 1 ) ) 

T 

74  . 

VEE 1 =VEE ( 1 ) *CONVl 

T 

74 . 1 

TIME  =  IT IME 

T 

75  . 

SL0RE  =  S1N (BOLONG-WF*T I  ME ) 

T 

76. 

CLORE=COS(BOlONG-WE*T I  ME )  _ 

T 

77, 

X  1  =C 1 *ClORE-C2* SLOPE 

T 

78. 

Y1 =C2*CL0RE+C 1 *SLORE 

T 

79. 

SLAR=C3*CvEE+C4*SvrE 

T 

80. 

CLAR=SORT ( 1 .-SLAR**2 ) 

T 

81  . 

c 

USE  SHEBO  TO  RESOLVE  PROGRADE  AND  RETROGRADE  MOTIONS 

IF ( SLAR ) 766.777 . 777 

T 

82. 

766 

IF (SHE BO) 56. 55 .55 

T 

83. 

777 

IF ( SHE BO) 55. 56 . 56 

T 

84  . 

55- 

T 

GOT 05 7 

T 

. 

86  . 

56 

C0NGR  =  ASN0DE  +  AC0S (  ( COM*CVEE-SOM»5VEE )/CLAR ) -WE*T I  ME 

T 

87. 

57 

CONT I NUE 

T 

88. 

ALT=P/( 1 .+ECC*CVEF) 

T 

89. 

Alt  1 =alt-rf 

T 

90. 

c 

REMOVE  periodicity  EFFECTS  FROM  VEE 

T 

ion. 

32 

1F( VEE 1-360. 0)30. 30. 31 

T 

101 . 

31 

VEEl=VEEI-360.0 

T 

102. 

GO  TO  32 

T 

103. 

30 

CONT 1 NUE 

T 

104. 

CLOR=COS(CONGR ) 

T 

105. 

SL0R=S1 N(C0NGR ) 

T 

106. 

CNUMP=  (  ITIME-IT  IM'^O) 

T 

107. 

RLAT= ASIN ( SLAR ) *CCNV 1 

T 

108. 

RLONG=CONGR*CONV I 

T 

I  09. 

JT I ME=T1MF/60. 

T 

109.1 

77 


Ti  me=t  iME'/aeon.n 

C _ PESOLVE  longitude  INTO  EAST  AND  WEST 

36  IF { RLONG ) 33» 34* 35 

33  IF { RLONG+ 1 00*0 ) 37, 34 t 34  _ 

37  RL0NG=RL0NG+360 *0 

GO  TO  36  _  _ 

35  IF { RLONG-I eoVo ) 34 , 34 ,38 


T  I  I  0* 


T  Ill, 

J-A  *  2* 

T  I  I  3  • 
1  14. 
T  115# 


38  RLONG=RLONG-360 *0 


T  1  16, 


c 

c 

GO  TO  36 

34  CONTINUE  _  _ 

REMAINING  STATEMENTS  CONTROL  RUNNING  OF  PROGRAM 
subroutines  for  integration  of  true  anomaly  and 

AND  CALLING 
calculat ion 

in 

OF _ 

T 

T 

1  I  7, 

I  1  0  •  _ 

c 

parameters  between  station  and  satellite 

NI N2sNUM{ 2 ) 

T 

I  18.  1 

NI N4=NUM { 4 ) 

IF{NUM{7)  )  70 1  #70 1  ,720  _  _ _  _  _ 

701  CALL  CALC  { NuM  •XI,YI,2I  • RE ,  1 T 1  ME , CL AR , SlAR , CLOR , SLOR , RL AT , 

2RL0NG«AlTI ,ALT, VEEJ ,BETA,SNODE)  _  _ 

I F { NUM ( 7 )  ) 72  I  , 72 1 , 720 
720  IF(NPRINT-H  )702t703,702 


I  18 
1  I  9 


120 

120 

120 


702  PRINT  61 

NPR1NT=-1  _ 

703  PRINT  65« JTiMEfTIME ,RLAT,RL0NG«ALTI ,VEEI 

721  GO  TO  (722, 7l4t  715  )  ,  1 _ _ 

722  I F { NUM { 6 ) )  724 , 724 , 723 

72  3  I F (  I T I ME-NUM ( 5 )  )  2 , 8 , 8 


120 

120 

120 

121 

122 

123 


724 
3 


714 

5 


IF{ I T I ME-NUM (3) )2,3,3 
CHECK=ALT 
I  =  I  + 1 
GO  TO  5 

IF { alt-check )6«6  *5 
DO  4  J= I ,N1N4 


124 

125 

126 

127 

128 

-129 


715 

2 


I 


CALL  RUNG  ( VEE { 1 ) , P , ECC , NUM { 0 ) ) 
4  CONTINUE 

I T I ME= 1 T I ME+NUM4 
GO  TO  700 

6  I  =  l-H _ 

GO  TO  2 

IF  { alt-re )0,8, 2 
DO  1  J= 1 ,NI N2 

CALL  RUNG  ( VEE { 1 ) , B , ECC , NUM ( 0 ) ) 
CONTINUE 
ITJME=JTI ME+NUM2 


130 

131 

132 

133 

134 


135 

136 

137 
136 

139 

140 


GO  TO  700  T  141 

44  FORMAT  {  I  1X,29H0RBI  TA.L  PARAMETERS  AT  BURNOUT/SX  ,  I  2HECCENTR  I  Ci  T  Y  a3X  i  T  141 
2 1  1 H 1 NCL 1  NAT  I  ON, 6X ,  1 2HSEM 1  LAT  REC • 4 X ,  1 2HTRUE  ANOMALY , 4X , 5H0MEGA ,  T  141 
36X,6HASN0DE/24X,3HDEG,  I_3X.f4HN  MI  ,  I  3X  ,  3HDEG  •  1  OX  •  3HDEG  •  9X  ,  3HDEG/7X  ,  T  141 
4F7.6,7X,F7.2,9X ,F0.2,9X , F 7 • 2 » 6X , F 7 • 2 , 5X , F0 . 3// )  T  I4l 

6 1  F0RMAT(9X,4HTIME^5X,4HTlMEt4X,69HVEHlCLE  LATITUDE  VEHICLE  LONGITT  141 


2UDE  VEHICLE  ALTITUDE  TRUE  AN0MALY/9X ,  I 4HM 1 NUTES  HOURS •  I  OX , 5hDE T 

3G  *  14X,aHDEG^6XtAtANt  Mi  j_12LKj  3ilQEG  1  _ T_141^7  _ 

65  FORMAT ( 8X ,  I  5# 2X ,F7 ,2 • 7X ,F8,3 »  1  1 X ,F8, 3, 1  OX ,F0,2 • 9X ,F8#3 )  T 

69  FORMAT ( IHI , 33X,50HSPACE  VEHICLE  TRACKING  AND  SENSOR  COVERAGE  PROGRT 
2  AM//  ^  -  -  -  ■  ■  -  -  ^ 

72  format ( //I  1 X,  IShBURNOUT  CONDI TI0NS/5X,8HALT1TuDE,3X,8hVEL0CITY,  T 


141 


141 

141 

141 

141 


9 

9^ 

90 

91 


26X,  5HGAMMA,4X»7HHEA0ING»4X,8HLATITUDE»3X,9HL0NGITUDE«4X,39H  SPIN  T  I  41 
3  AX  I  S  ATT  I  TUQE  SPIN  AX  IS  HE  AD  1  NC/7X  ,  4HN  Ml  ,  6X  ♦  6HFT/SEC  !»  PX,  3HDEG  •  T  14  1 
47X,3HDEGf8X,3HDEG,9X,3HDEGf I5X,3HDEG, 1 8X , 3HDEG/5X ,Fe • 2 • 3X ,F8, 2 ,  T  141 
54X,F7.2,3X,F7,2  t 4X , F7 . 2 t 5X , F7 • 2 » 9X , F9 ♦ 2  t 1 4X,F7,2//)  T  141 


92 

93 

94 

95 


78 


87  format (  I  IX t 25HPP0GPAM  CONTROL  CONST ANTS/5X «  1  1 0HNUM ( 1  )  NUM(2)  NUM ( 3T  1 4 1  • Q6 
21  Ny^^(4)  NUM(F)  NUM(6)  NUM(7)  NUM  (  o  )  KrUM(9)  NUM  (  1 P  >  NUM  (  1  1  )  NUM(12T  141>97 

3)  NUM(13)  NUM(I4)  NUM ( I  5 )  NUM (  1 6  ) /6X ,  I  3 » 4X ,  I  3 , 2X ♦  I  7 , 2x •  I  3 • 2X •  I  7 • 3XT  14  1.90 

4.!2t4X.  130X<14.4X<  I3.4Xtl3.6(5X.  13)  )  _  _  _  _  _T  J41.99 

8  RETURN  T  142. 

END  _  _  _  _  _ T 
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> 


*. 


! 


1 


APPENDIX  III 


FORTRAN  IV  LISTING 
SUBROUTINE  CALC 


c 

SYMBOL 

UNITS 

DEF IN  I T 1  ON 

c 

alat 

(DEG) 

LATITUDE  OF  ITh  GROUND  SENSOR 

c 

ALONG 

(DEG) 

longitude  of  ITh  GROUND  SENSOR 

c 

ALPHA 

(DEG) 

angle  between  line  of  Sight  and  ith  fixed  ground 

c 

SENSOR  ANTENNA 

c 

anaz{ 1 : 

)  ( DEG ) 

azimuth  of  ith  ground  sensor  antenna. if  fixed 

c 

;'>NEL(  I  ; 

I  (DEG) 

elevation  of  ith  ground  sensor  antenna, if  fixed 

c 

NGLE 

(DEG) 

ANGL^  BETWEEN  L I NE  OF  SIGHT  AnD  SPIN  STABILIZED 

c 

SPACE  BORNE  ANTENNA 

•  c 

AZD 

(DEG) 

AZIMUTH  ANGLE  OF  LINE  OF  SIGHT  BETWEEN  I TH 

c 

GROUND  SENSOR  AND  SPACE  BORNE  VEHICLE 

c 

BW  1 

(DEG) 

BEAMWIDTH  OF  THE  SPACEBORNE  AnTENNA  MEASURED 

BW( n  tQEGj 
DBA { I )  (NONE } 
ELE  (DEG) 


BETWEEN  HALF  POWER  POINTS 

beamwidth  of,  iihlorqund. S£nsqr  measured  between 
HALF  POWER  POINTS 

_DUMMY  INPUT  variable » set  equal  TQ  ON-AX I S  GAIN  OF 

the  I th  ground  station  tracking  antenna 
elevation  angle,  of  line  of  sight  between  ITH 


ELEM(  I  )  (DEG) 


FREO 

FNGLE 


5  ) 

(DEG  ) 


GROUND  SENSOR  AND  SPACEBORNE  VEHICLE 
MINIMUM  ELEVATION  ANGLE,  OF  I Th  GROUND  STATION 
TRACKING  ANTENNA 
OPERATING  frequency 

angle  between  line  of  sight  And  earth  pointing 
Space  borne  antenna 


ITIMEO 
RhO 
XAN (  I  ) 
YAN{  I  ) 
ZAN (  I  ) 
XF 


(SEC)  time  at  preceeding  integration  interval 
(N  MI)  _ line  of  _SIGHT  DISTANCE  BETWEEN  ITH  GROUND  SENSOR 
(NONE)  DIRECTION  COSINES  OF  ITH  GROUND  SENSOR  WITH  FIXED 
_  ANTENNA  AXIS 

(NONE)  directional  COSINES  FOR  SPACE  PORNE  ANTENNA 


C  YF  POINTING  TOWARD  EARTH  CENTER 

C  7F  _  _  _ 

C  THIS  SUBROUTINE  USES  CALC  TO  MAKE  IT  COMPATIBLE  WITH  OTHER 

c  calling  progams  _ 

Subroutine  calc  (num,x i • y i • z i  •re* i t ime*clar,slar»clor,slor»rlat,c  i* 

2RL0NG*ALTI •alt* VEE 1 •BETA^SNODE)  C  1*1 


DIMENSION  ALAT (150)  • ALONG (  I  50 )  •CL A (  150)*  SLA (  I 5O  )  •CLO (150)  • SLO (  I  50 )C  1.2 
2 •ROH{  150 )  •EEL (  I  50 )  • ANGEL (  150 )  •ADZ (150)  •FNGEli  I  50 )  • AnEL (150)  • ANAZ (  I C  1.3 
350 )  • XAN (  150)  • YAN(  I  50 )  • ZAN (  150)  • BW (  150)  • NUM (  I  6 ) • DB A (  I  5 0 )  • ALP AH (  1 50 ) C  1 

4tALTS(  150)  •ELEM(  150)  •SNAmE(2*  150) _ _ _  C  1.401 

C0NV=6. 01 745329  '  C  l.4l 

CONV  I  =57.29576 ^ C  I  .42 


CFEET=6076«  I  ! 

T I  ME- I T 1  ME 

JT I ME  =  T I  ME/60.0 

! I ME=T IME/3600. 0 

READ  IN  GROUND  AND  VEHICLE  SENSOR  DATA 
IF ( MSTART+I )II6*II5*1I6 


I  .43 
1  .44 
I  .45 
1  .46 


1 1 6  MST ARTS- I . 

PRINT  73 

READ  195*FRE0*BWI 
195  FORMAT (2F 1 0.2 ) 


1  .6 
I  .61 
1  .7 
1.71 
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BWG=BW1 

c 

1  ^7P 

BW  1 =BW I ♦CONV 

c 

1 .73 

NUGS=NUM ( 1 ) 

c 

1  *8 

DO  lOl = 1 .NUGS 

c 

2. 

READ  78 .  I DUM ♦  ( SNAME (J.I)«J=I.2)«  ALAT (  I  )  , ALONG (  I 

) . ANEU ( 1 ) 

. ANAZ  f 

3. 

2ALTS( I ) .BW( I ) .ElEMI I ) 

3.  I 

78  FORMAT  ( 2X •  I  3 • 3X • 2 A8  .7F8.2) 

C 

3.2 

PR  I  NT 75  .  I  •  ( SNAME ( J.  I  )  . J= I  <2 )  •ALAT (  I  )  . ALONG (  I 

)  ♦ ANEL (  I  ) 

•  ANAZ (  1  )  •  C 

3.2 

2AlTS( 1 ) .BW ( 1 ) ,ELEM( I ) 

C 

3.25 

75  FORMAT  ( 2X .  I  3 • 3X . 2 A8  .F7.2.2F9,2 

• IX.F9.2. 

3X.F9.2. 1 XC 

3.3 

2 .F9 .2 .9X .FT. 2 ) 

C 

3.4  . 

CLAS=CONV#ALAT ( I ) 

C 

4.3 

CLOS=CONV*ALONG ( 1 ) 

c 

4  .4 

CLA ( I ) =COS(CLAS ) 

C 

4.5 

SLA( I )=SIN(CLAS) 

c  . 

CLO(  I  )*COS(ClOS) 

c 

4*7 

SLO( 1 ) *SIN( CLOS ) 

_  c 

4.8 

ANELR=CONV*ANEL ( I ) 

c 

4.9 

ANAZR=CONV*ANAZ ( I ) 

c 

5.0 

CANEL=COS( ANElR ) 

c 

5.2 

SANEL=SIN ( ANELR ) 

c 

5.3 

CANAZ=COS( ANAZR ) 

c 

5.4 

SANAZ=SIN( ANAZR ) 

c 

5.5 

CAl =SANEL*CLA( I ) -CANEL*CANAZ*SLA ( I ) 

c 

5.6 

CA2  =  CANEL»SANAZ 

_ 

^LtJ 

-±D. 


XAN(  I  ) =CAI *CLO(  I  ) -CA2*SLO(  I  ) 

ZAN( I )=SANEL*SLA( I ) +CANEL»C ANAZ*CLA ( I ) 
XQN-LINUE 


I  15 


XF*CLAR*CLOR 

YF=CLAR*5LQR 


C 

C 

C 

C 

C- 


5.8 

5.9 

6.0 

- 

6.11 

6.1Z 


zf=slar 

M=0 


c 

x 


6.  13 


DO  62  L* 1 •NUGS 
Ai  T,c;isAl  T5^fL>yrFFrT 


C 

-X- 


7. 

7>  1 


CPSI »SLA(L>*SLAR+CLA(L)*CLAR*(CLO(L> ♦CLOR+SLO(L) ♦SLOP) 
RHQx^QRT(ALT#»2-KRF+ALTSM»»2-a^»ALT»tRE-^ALTS!  1»CPSI  \ 


C 

C 


8. 


SELE* ( ALT^CPSI- (RE+ALTS I ) ) /RHO 

CELEgSQRTM  .  -SELE*»2 )     

ELE»ASIN(SELE) 

-EtLL  =(SL^R-SLA(H»CPS1  )  /  (  C1_A  (  L)  »SQRT  (  1  ,  fJtZXA- 

AZR= ACOS( PHI ) 

AZD  =  A2R»CQNV1 


C  I 

C  I 

C 1 


C  1 
C  1 


4, 

ix 


C 

X- 


SWITCHING  TO  RESOLVE  AMBIGUITY  FOR  EAST-WEST  LAUNCHINGS 

And  south  latitude  mode _ _ _ _ 

clongr*rlong 

C1-QNGA  =  ALQNG(L) 


IN  NORTH 


C  16. 

X  -IX.—  - 


IF (CLONGR-CLONGA)88,92,90 
88  1F(C1 QNGR-C1  nNGA^>  1  80  >  O  \  9P  ,  Q?  .  Q7 


C  18. 

c  la^. 


90  IF (CLONGR-CLONGA-I 80.0 )92*92*93 
93  AZD=360.0-AZn 
AZR= AZD*C0NV 

92  CON  T  I  NUE _ _  _  _ _ 

CAZf>*GO&4  AZR» 


C  20. 
-CXI.. 
C  22. 
C  X3  .-_ 
C  -24  1- 


SAZR=SIN( AZR ) 

C5^SE4X*CLAtL  )  -C£UE  »CAZR»SLA  (  L  ) 
C6=CELE*SAZR 

XRHOiC5tCLQ<U^-C61^SLO<U> 
YRH0*C6*CL0(L)+C5*SL0(L ) 

-ZRHQ  =  SELF#S1  A  f  1  1 ->-CFl F  #GA7P#  n  Aft 


25. 

26  . 

27. 

28. 

29. 

30. 


ANGLE^ACOSC  X  I *XRH0+Y I »YRH0+Z I »ZRHO) 
FNGLFsArQS<yF»YPHn4-VF#YRHQ^-ZF»7PHn)  
alpha = ACOS ( X AN ( L ) ♦XRHO+YAN (L ) »YRHO+ZAN( L ) ♦ZRHO ) 
--  Fi  F=ELEtCQNVl  _  „  .  _ 


C  31  • 

c  aa. 

C  33. 
C  34. 
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ANGLE=ANGLE*CONV I  C  35* 

_ FNGL  E=FNGLE*CQNVJ _ C  _ 

AlPHA=ALPHA*CONVI  C  37* 

C  GAl  rULATF  JlVEgAGE  RATEL- OF CHANGE!  OVER  CAST  PRINTOUT  INTEiiVAL- 

CNUM8= I T I  me- I T I MEO  C  38 • 

_ lEL.(l.STARliLJJ_151jl5Qj-i51 _  _ _  _  _ 

150  DRH0= ( RHO-ROH ( L ) ) /CNUMB*CFEET  C  39* 

nAznr  ( A7n^An7  (L  1  i/r,NUMR - C  .tiQ^ 

DELE=  (ElE-EEL  (  L  )  ) /CNUMB  C  i  . 

nANG<  -F  >  (  ANGl  f^anGEJ_(LA^7.CNUM.3 _  _  _  _ _  _  C  42. 

DFNGLE= (FNGLE-FNGEL (L) ) /CNUM8  C  43* 

_DALPHA=1ACPHA=ALPAH1J_T  )V-CNUMa  _  C  44. 

151  LSTART=-1  C  44*1 

^ _ SAVE  VALUE  OF  VARIABLES  AT  PRECEEDlNG  mTHRVAi _ 

ALPAH ( L) =ALPHA  C  45* 

_ROH(Lj=RHO _  _  _ _  ,  C  46. 

EEL(L)=FLF  C  47. 

ANGEL  (L  )  =  ANGLE  _  C  48, 

FNGEL(L  )=FNGLF  C  49. 

ADZ(L)  ^AZD_ C  .50^ 

K=L  C  51 . 

C  _  EL  1  MIN  A  TE  DAT  A  WITH  ELEVATIONS  BELOW  ACCEP  T  ABlE  M  I  N 1  M.U  M 

IF ( ELE-ELEM ( L ) ) 50«64 *64  C  51.1 

50  M  =  M+J  _  _  _  _  C  51*2 

iriNUCS  M;b:«5!.62  C  L1.3 


_ 51  ALPHA  =  0. _ C  51.3  1 

ANGLE=0.  C  51.32 

_ FNGLE  =  0. _  _ ^  C  51.30 

IF( JSTART+I )64«30«64 

64  1F(NUM(9)  )  1  1  1  <  1  1  1  <  1  12 _ _  _  C  52. 

C  CALL  subroutine  DEBE  TO  CALCULATE  LOSS  DATA  IN  DB 

111  CALL  DEBE  (  BW  .  ANGLE  .  ALPHA  .FNGLF  ■ _ BW  1..BHO.K  .DBA  »  ,SLOS  . _ C  53. 

2D81  •DB2«D83fD84  ^ DBS •DB6«DB7«D08  fBWG* DBG « FREQ* GANGLE«GALPHA,GFNGLt ) C  54 . 

,  112  C0NI_1NUE!_ _ _  _  _  C  5A^  i 

IF( JSTART+1 )lAlfl40«l4I 

_ 14  1  ^STARXS^J _  _  _  -  c  ^.4 

PRINT  76«BWGfDBG*FRE0  C  56. 

_ PRlNXJSfi _ _ 

,  140  IF ( KSTART+ 1 ) 143 • 1 42* 1 43 

o _ 143  KSTART=:-1 _  _  „  q  57^1 

•  30  PRINT  69* JT IME* T IME •RLAT fPLONG* alt  I  • VEE 1  •BETAfSNODE  C  58. 

*  IF(NUGS-M)_6ai62, 144  C  52^01 


144  PRINT  BF 


-C- 


PRINT  89 

142  PRINT  60*  (SNAME( J*L) • J=1 *2) • ELE • A2D • RHO • ANGLE • FNGLE « ALPHA , GANGLEC 
2*GALPHA.DB1  • DB3 • DP5 • DB7 • DELE • OAZD • DRhO • DANGLE • DFNGLE  <  DAL^H A . GFNGLFC 
3 • SLO S • D82  » D84 , DB6 , DBS  C 

62  CONTINUE  C 

KSTART  =  >1 _ C_ 


58^4 

58.2 

59. 

60. 
61  . 
62. 


60  FORMAT  < 2X • 2 A8 • 3X • F 7 . 2 • 2X • F7 . 2 .  1 X • F9 . 2 • 4X • F8 • 3 • 4X • F7 . 2 • 2X • F7 . 2 • 5X •  C 
2F7^2v^)^jF7.2j  IX  ,E7.2l22IXtF7^^2X.F7^4  1 1 X. F9 .  1  «  4X , F8 •  4  •  4X,F7.4.  C 
32X*F7.4,5X*F7.2 •2x •F7.2  ^4 (  1X«F7.2 )/)  C 

76  FORMaT<//10X.29HVEHICCE!^ANQ^SENSQR  PASAJ^ET£RS/5X «  10H8E AM  WIDTH, 3X,C 
28H TRACK ING. 5X • 9H0PERAT I NG/8X  «  3HDEG • 9X , 4HGA I N, 7X , 9HFREC 

_ 30yENj;_y _ _ /5X^E7^2.5x.f7.p  .  tv  .fa.?  i  r 


64. 

65. 

65.  1 

66. 
67. 

— 


73  FORMAT! 1 lX,aH5TATlONS/lX,6HNUMBER, IX,4HNAME* 15X,3HLAT,5X,4HL0NG«  C  70. 
25X  ,  4HEL£y,*_5X  ,  ^HAZgMUTH  t  3X  tOKACTi  TUOBjtBX^OHBEAM  W  IDTH , 3X  .^1 4HM.1N  ELC  71. 
3EV  ANGLE/27X,3HDEG,6X,3HDEG«6X,3HDEG,7X,3HDEG,7X«4HFEET,8X<  C  71.1 

4  3HDEG  ,  I  3X  ,  3HDEG  ) _  C  71  .2 

86  FORMAT < //I  X. 62HT0TAL  LOSSES  FOR  VARIOUS  VEH I ClE-GROUND  ANTENNA  COnC  72. 
2F1GURAT  IONS  C  73. 
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2  /5Xf46HDBlsVEHlCLE  SPIN  STABILIZED-GROUND  FIXED  ♦  C  74, 

2  /5Xf  46HDB2  =  V£HICLE  SPIN  STABILIZED-GROUND  TRACKING _ i _  C  75  > 

2  /5X*46HDB3=VEH1CLE  EARTH  CENTERED  -GROUND  FIXED  »  C  76. 

_  2  /5X*46HDB4=VEHICLE  EARTH  CENTERED  -GROUND  TRACKING  t  C  77t 

2  /5Xf 46HDB5«VEH ICLE  TRACKING-GROUND  FIXED  » 

2  /5X  <46HDB6«VEHI  CLE  TRACKING-GROUND  TRACKING _ j _  !  ^  ^ 

2  /5Xt46HDB7  =  VEHlCLE  SPIN  ST AB 1 L 1  ZED- 1 SOTROP I C  GROUND  »  C  9C . 

2  /5Xf 46HDB8gVEHlCLE  EARTH  CENTERED- I SOTRQP 1 C  GROUND  )  C  81. 

69  FORMAT!  // 1 X ♦ 7H VEH I CLE f 4X . 4HT I  ME . 6X . 4HT I  ME t 4X . 7 1 HLAT I TUDF  LONG  1 C  82. 

2TUDE  ALTITUDE  TRUE  ANOMALY _ OMEGA  _  A  SNODE / 1 X  v  8HP0^1  T  C  83 . 

310Nf 2Xf7HMlNUTESf4Xf5HHOURSf 5Xf3HDEG. 10Xf3HDEG» 10X.2HNM. 1 1X.3HDEG.C  84. 

4  10Xf3HDEGf 1  OX . 3HOEG/ 1 2X .  1 5  * 4X . F6 . 2 . 2X . F8 . 3 . 5X . F8 . 3 . 5X . F8 . 2 f 3 < 5X . F8C  85. 

5.3))  C  85 . 1 

85  format (/2XfeHSTAT 1 ONS f 9X . 9HELEV AT 10N.2Xf 7HAZ I MUTh . 5X . 5HR AN GE * 3X . 9 HQ _ 86* 

2 ant  angle f2X »9H ant  ANGLE . 2X . 7HGRD  ANT . 1 X . 1 1 HVEH 1 CLE  ANT.2X.7HGRD  AC  87. 
3NT.2X.29HVEHICLE-GROUND  ANTENNA  LOSSES/2X t 2H I N ♦ 1 5X t 9HELE V-RATE f 2X . C  88. 
47HAZ-RATE  *1X»9HRANG  RATE • 3X . 9HSp 1 N-ST AB f 2X f 9HVERT-ST AB » 4X . 5HANGLEC  88.1 


5»4X.8HGA1NS  DB*2X.7HGAIN  De/2X < 7HCONT ACT . 1 6X . 3HDEG . 6X » 3HDEG . 3X .  C  68.2 

67HNAUT  M 1 .RX f3HDEG»8Xf 3HDEGf 6Xf 3HDEG. 3Xf 9HSPIn-STAB»4X.5HSPACE ♦  C  88.3 

75X.3HDB1 f 5X.3HDB3»5Xf 3HDP5*5Xf 3HDB7  )  _  _  _  C88.4 

89  FORMAT (21X.7HDEG/SECf2X.7HDEG/SEC.2X.8HFEET/SEC»5Xf7HDEG/SEC«4X.  C  88.5 

2  7HDEG/SEC  *2X  f  7HDEG/SEC<  3X  .  9HVERT-ST  AB  *  2X  .  7HLOSS  DB  »  5X  .  3HDe2^f  _  C  88.6 

35Xf 3HDB4f5Xf 3HDB6.5X.3HDB8/)  C  88.7 

ITIMEO-ITIME  C  94. 

RETURN  C  95. 

END  _ _ _  _  _  C  96. 
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APPENDIX  IV 


FORTRAN  IV  LISTING 
SUBROUTINE  RUNG 


^ _ This  is  SUtJRQUTlNF  RUNG! _ 

C  THIS  subroutine  USES  RUNG  TO  MAKE  IT  COMPATIBLE  WITH  OTHER 
C  calling  PROGAMS 

C 


c 

subroutine  rung  (VEE  iB*ECCiNUMB> 

DIMENSION  C<4 ) ♦ VEE ( 5 > . VEER ( 4 ) 

A  MODIFIED  FOURTH  ORDER  RUNGE  KUTTA  NUMERICAL 

INTEGRATION  METHOD 

R 

R 

I  . 

2. 

c 

IS  USED  IN  THIS  SUBROUTINE  TO  OBTAIN  THE  TRUE 

anomaly  as  a 

c 

FUNCTION  OF  TIME  IN  SECONDS.  NUMB  IS  INTEGRATION  INTERVAL  IN  SECS 

. 

IF  <  I  START + I  >104,103.104 

I  04 

C<  I  >  =0.5 

R 

4  • 

C<2)=0.5 

R 

5. 

C(3)  =  I  .0 

R 

6  . 

C<  4 ) =0.0 

R 

7. 

TNUM=NUMB 

R 

7. 1 

ISTART=-1 

R 

6  . 

103 

D0I02K=I ,4  _ 

R 

9. 

VEER<K  )  =B*(  I  .+ECC*COS ( VEE <K )  )  ) ♦♦a^TNUM 

R 

1  0. 

I  02 

VEE(K+I  )=VEE(  1  )-^VEER(K)*C(K) 

R 

1  1  . 

I  0  I 

VEE  (  I  )  =VEE(  1  )  +  (  VEER  (  1  )  +  (  VEER  (  2  )  -♦-VEER  (  3  )  )*2.-^-VeER  (  4  )  )/6.0 

R 

I  2. 

RETURN  _ 

R 

13. 

END 

R 

14. 
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APPENDIX  V 


FORTRAN  IV  LISTING 
SUBROUTINE  DEBE 


SUBTY.be  4FOQICAM*LMAg 


C 

C 


list  of  symbols  used  in  subroutine  DFMEl 


c 

SYMBOL 

UNITS 

DEF INI T ION 

D 

G 

Angle 

(OEG) 

angle  between  spin-stabilized  SPACEBORNE  ANTENNA 

0 

C 

AND  THE  LINE  OF  SIGHT  TO  THE  ITH  GROUND  STATION 

D 

Q 

Alp*--*  a 

^  np-r-.  ^ 

Atwini  p- FlgTujFg-Kj  Tt-IF- Av  I  C  QT  Tur  F  I  vPn  AKITFMKiA  OF  TMFO 

c 

ITH  ground  station  AND  THE  VEHICLE  LINE  OF  SIGHT’ 

D 

c 

BW 

BFAMWIDTW  of  the  ITH  GROUND  ANTFNNA  MEASURED 

D 

c 

BETWEEN  THE  HALF  POWER  POINTS 

D 

G 

BwG 

(DEG) 

Dummy  variable  used  to  preserve  bwi 

O 

c 

Bw  1 

(DEG) 

BEAMwIDTH  OF  THE  SPACEBORNE  AnTENNA  MEASURED 

D 

r 

RFTbjpFN  WAl  F  DOWFD  OOImTS.   .  

n 

C 

Q 

DBA 

(  DB  ) 

ON-AXiS  gain  of  THE  ITH  GROUND  STATION  TRACKING 

D 

o 

c 

DBG 

(DB) 

On-axis  gain  of  the  space-horne  tracking  antenna 

D 

iC 

DSl 

XOBI 

OF  THE  SPIN  ST Ati  1 H ZED-GROUND  FIXED 

D 

c 

ANTENNA  CONFIGURATION 

D 

r 

nnp 

( nR> 

TnTAl  1  nSS  DF  THF  SPIN!  STAHIl  tZm-GROUND  TRACKING 

n 

C 

antenna  CONFIGURATION 

D 

c 

DB3 

toa_i 

TOTAU  LOSS  DF  THE  EARTH  CENTERED -GROUND  FIXED 

D 

c 

ANTFNNA  CONFIGURATION 

D 

cu 

DBA 

(DB) 

TOTAL  JLQSS -OF  THE  £J\RTH  CENTERED-GROUnD  TRACKING 

D 

c 

ANTENNA  CONFIGURATION 

D 

c 

DBS 

(DB) 

TOTAL  LOSS  OF  THE  VEHICLE  TR ArK I NG-GROUND  FIXED 

D 

c 

antenna  configuration 

D 

.PB6  <QB)  total  LOSS  PF  TH£:  VEHICLE  TR ACK I NG“GR0UNQ  TRACKINGD 

ANTENNA  CONFIGURATION  D 

_DB7  .  _TQT,AL-L.PSS  QF_  the:  SlP IN  STAB  1 L  I  ZED- 1  5ENTROP 1 C  D 

GROUND  ANTENNA  CONFIGURATION  D 

DBS (NONFl  TOTAL  LOSS  OF  THF  FARTh  CFNTFPFH- 1  SFNTPOP  t  r  f^rVv^LiNnn 


C 

JL. 


c 

PDSL 

(D6) 

C 

C 

FNGLE 

(DEG) 

C 

FREO 

(MCS) 

c 

GALPHA 

(DB) 

C 

c 

c 


ANTENNA  CONFIGURATION 

FREQUENCY  DEPENDENT  _PART  OF  SPACE  LOSS  D 

angle  between  earth  center  POINTING  SPACEPORNF  D 
antenna  and  line  Q£_S1GMT  of  ITJi  G.ROUND  STATION  D 
OPERATING  frequency 

GAIN  OF  _THE  F I XED  ANTENNA  OF  ITH  GROUND  STATION _ D_ 


C 

-C 


GANGLE 

GFNGLE 

L 

RDSL _ 

RHO 


(DB) 

(DBl 


( NONE ) 

(DBJ 

(  N  MI) 


GAIN  OF  THE  SPIN  STABILIZED  SPACE-BORNE  ANTENNA 
-GA,iN_oe_IHE .EARTH  centered  SPACErBORNE  ANTENNA 
STATION  NUMBER  OR  INDEX 
J?ANGE  DEPENDENT .  PART  OF  SPACE  EOSE 
LINE  OF  SIGHT  DISTANCE  BETWEEN  THE  ITH  GROUND 
STATION  AND  THF  c^PArF-RORNE  VFHIGLF 


D 

D 

D 

D 

D 

JZL, 


C 


SLOS 


(DB) 


TOTAL  SPACE  LOSS  (FDSL+RDSL) 


C  THIS  subroutine  CALCULATES  THE  SPACE  LOSSES  AND  THE  GAINS  OF 

_C _ various  vehicle-ground  AT  ATI  on  ANTENNA  CQNFIGURAT I  QNS*  ALU  INPUTS 

C  ARE  PASSED  THROUGH  THE  CALL  STATEMENT  IN  SUBROUTINE  CALC# 

C  THIS  SUBRQUX1nE USES  DEBE  TO  mAkF  IT  CDmpatiHlp  WITH  hthfp 


c 

c 


calling  progams 

subroutine  DEBE  ( BW • ANGLE • ALPHA • FNGLE •  BW 1 • RHO • L • DB A • SlOS *  D 

2DB I  • DB2 . DB3^Q04 • PB5*DB6 • 0E7  • OBQ . BWG  *060 . FREQ • GANGLE  t  GALPHA i GFNGLL ) D 
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PiM£N6iW  iSO  >  ^OHA<  150) - 

CONV  =0.01745329 

CONI  a  c3,54491/  caW-(UMCONV  ).)*Jt2*0.6 
C0N2=(3.‘=>4  491/BW1  )  »»2#0.8 

FDSL=2D.*AL0G1QJ  (Z99._776*3.ZflOa3^6076,  1  15)/FREQ-) 
RDSL»-20. »AL0G1 0 (4 •♦3. 1 4 15927*RH0 ) 
PSG=1Q.^>AL0G1Q(C0n2) 


-«?-• — 

3. 

3.5 

*^.6 

.  1 

4  .2 

-.a-. 


o 

D 


(' 

JX- 


124 


125 


DBA ( L)  =  1 0 • * ALOG 1 0 ( CON 1  ) 

IF  (ALPHA) 125. 124^ 125  _ 

ARC  =1.0 


GQ  TQ. 

ANl 

ARG 


IZG _  -  _ 

=2. 78332* ALPHA /BW(L) 

s  <  g^IMt  AMI  \yAKii 


II. -It*?  ■ 


D 

O 

o 

a 

o 

-0- 


1 1  • 

13. 

14. 

15. 

16. 

O  18. 
£y  I9w 
D  20* 

a  21^ 

D  22. 

O  23. 


126  GALPHA= 10.*ALOG1 O(C0Nl*ARG) 
_ LF  (ANGLE  1135.  134.135-  -- 


1  34 


ARG2 
GO  TO 


=  1  .0 

136-- 


135 


AN2 

ARG2 


=2.78332*ANGLE/BWG 

=  <c;iMf  AMP>yAMP>*#P 


1  36  GANGLE=  1  O'. ♦ALOG  1  0  (  CON2*ARG2  > 

234  ARG3  =1.0 

_ GQ  .TG  ^36 _ _ _ 

235  An3  =2.78332*FNGLE/BWG 
ARG 3 =(S1N(AN31/AN3>**2 


24, 

25* 

26. 

27*. 

28. 


n  29. 


236  GFNGLE=10.*AL0G10(C0N2*ARG3) 
__  ADD  UP  LOSSES  TN 

SLOS  =fosl+rdsl 
_ DBI, 

DB2  =SLOS  +  GANGLE  +  OBA ( L  > 

DB3  «SLOS-^GFNGLE^^GALPHA 


SYGTE^i. ANTENNA  CONFIGURATIONS 


O  30. 

O  31  • 
D  32*, 
O  33. 
D  34. 


OB  4 
J)^5_ 
DB6 
OB  7 


=SLOS+GFNGLE+OBA ( L ) 
=  SLOS->-OBG^-G  ALPHA 

=  SLOS+DBG  +  DBA (L  > 

OS +G  ANGLE 


DB8  =SLOS+GFNGLE 
return 


O  35. 

D_ 

D  37. 

0  38.„_ 
O  38. 1 
D  39. 


END 


O  40. 
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APPENDIX  VI 


FORTRAN  IV  LISTING 
SUBROUTINES  ASIN  AND  ACOS 


T  subtype*fortran*lmap 

CASIN 

_ FUNCTION  _ _  _ _ _ 

DATA  RAO  (57.295779513) 

IF  (ABS(A)-I.)  I*2«2 

2  IF  (A)  20«2l*2l 

20  IF  ( A-f  I  .OOOOOI  )  22*23*23 

23  ^SlN*-90«/RAD 
GQ.  TQ-IO 

22  IF  (A+I.OOI)  24*25*25 
25  ASIN--90./RAD 
PRINT  30 

30  FORMAT(6aH  ARG  OF  ASIN  S I GN IF  I CANTL Y  GT  1.  ARC  SET  *  X  AND  PROCR 
2AM  CONTINUED) 

GO  TO  10 

24  ASlN»-90«/RAD 
PRINT  31 

31  FORMAT ( 6 IH  ARG  OF  ASIN  MUCH  GREATER  THAN  ONE  -  PROBALE  PROGRAMING 

2ERR0R) 

GO  TO  10 

21  IF  (A-l.QOQQQl)  43.43^42 

43  ASIN«90./RAO 
GO  TO  10 

42  IF  (A-l.OOI)  45*45*44 
45  ASlN«90t/RA0 
PRINT  30 

_ GQ  TO  1 Q _ 

44  ASIN«90./RAO 
PRINT  31 

GO  TO  10 

I  ASIN*ATAN( A/S0RT( I f 0“A*»2) ) 

I  0  RETURN 
FND 


T  subtype*fortran*lmap 

GACOS  - -  —  - 

function  ACOS  (A) 

DATA  RAO  (57.295779513) 

IF  (ABS(A)“-l.)  1*2*2 

2  IF  (A)  20*21*21 

20  IF  ( A-*-!  .OOOOOl  )  22,23*23 

_ ^3L^ACOS=19Qf/RAD _ _ _ 

GO  TO  10 

22  IF  (A+I.OOl)  24*25,25 
25  ACOS« I80»/RAO 
PRINT  30 

30  format  ;66ri  AkG  OF  ACOS  SIGNIFICANTLY  GT  1.  ARG  SET  »  I  AND  PROGRA 
2M  CONTINUED) 

GO  TO  lO 

1'4  ACOS=  ISC./RAO 
PRINT  31 

31  FORMAT (62H  ARG  OF  ACOS  MUCH  GREATER  THAN  ONE  -  PROBABLY  PROGRAMING 
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2  ERROR) 

GO  TO  10 _ 

21  IF  (A-1. 000001)  43«43f42 

43  ACOS-0.0 
GO  TO  10 

42  IF  (A-l.OOl)  45«45,44 
45  ACOS-0.0 

PRINT  30 _ 

GO  TO  10 

44  AC0S«0«0 
PRINT  31 
GO  TO  10 

1  1F<A)6»5*7 

6  ACOS»ATAN(SQPT< 1 • 0-A»»2 ) /A ) ♦ 1 60 »0/RAO 
GO  TO  10 

7  ACOS«ATAN(SORT( 1#0-A**2)/A) 

GO  TO  10 

5  AC0S-90.0/RAD 
10  RETURN 
END 
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APPENDIX  VII 


DERIVATION  OF  THE  ANTENNA  ASPECT  ANGLES 

To  determine  the  space  losses  and  antenna  gains  for  different  orientations 
of  ground  station  and  vehicle  antenna  axes,  it  is  necessary  to  derive  expressions 
for  the  angles  between  the  antenna  axes  and  selected  space  directions.  In  the 
DEBE  subroutine,  combinations  of  the  following  four  distinct  situations  are 
considered. 

1)  Ground- station  antenna  fixed 

2)  Ground- station  antenna  tracking 

3)  Vehicle -antenna  fixed 

(a)  Spin  stabilized 

(b)  Stabilized  on  earth  center 

4)  Vehicle-antenna  tracking 

The  angle  relations  between  selected  space  directions  are  simply  obtained 
by  resolving  unit  vectors  along  these  directions  into  components  along  the  X,  Y, 
and  Z  axes  of  the  reference  earth-fixed  coordinate  system  and  by  taking  the  dot 
product  of  the  appropriate  transformed  unit  vectors. 

RESOLUTION  OF  FIXED  GROUND-STATION  ANTENNA 

The  direction  of  a  fixed  ground-station  antenna  is  given  by  the  azimuth 
(ANAZ)  and  elevation  (ANEL)  angle  of  the  antenna  axis  relative  to  the  local 
horizontal  plane  (see  Fig.  15).  Clockwise  directions  from  north  define  positive 
azimulh  angles  and  the  up-direction  define  positive  elevation  angles.  Two 
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SA 


Fig.  15  Geometry  for  Unit-Vector  Resolutions 


resolutions  transform  a  unit  vector  in  this  direction  into  components  in  earth- 
centered  coordinates.  The  first  resolves  the  vector  into  components  along 
local  tangent-plane  axes:  X'  pointing  north;  Y'  pointing  west;  and  Z'  pointing 
along  the  radius  vector  from  the  earth  center.  The  second  resolves  these 
intermediate  components  into  the  components  XAN,  YAN,  ZAN  along  the  earth- 
centered  axes.  The  orientation  of  the  local  tangent  plane  relative  to  the 
earth-centered  axis  system  is  defined  by  the  geometric  latitude  (ALAT)  and 
longitude  (ALONG)  of  the  station.  The  geometry  between  earth  and  antenna  axis 
is  shown  in  Fig.  15.  The  required  resolutions  are: 
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From  the  unit  vector  e^^^  to  components  X'  ,  Y'  ,  Z'  along  local  horizontal 

axes 


X' 

cos  (ANEL)  cos  (ANAZ) 

®AN 

Y' 

= 

cos  (ANEL)  sin  (ANAZ) 

0 

Z' 

sin  (ANEL) 

0 

From  the  local  horizontal  axes  components  X'  ,  Y'  ,  Z'  to  earth-centered 
components  XAN,  YAN,  ZAN. 


XAN 

-sin(ALAT)  cos(ALONG)  -sin(ALONG)  +  cos(ALAT)  cos(ALONG) 

X' 

YAN 

-sin(ALAT)  sin(ALONG)  +cos(ALONG)  +  cos(ALAT)  sin(ALONG) 

Y' 

ZAN 

cos  (ALAT)  0  sin(ALAT) 

Z' 

(62) 


Combining  these  matrices  gives  the  following  components  of  the  unit  vector 
along  the  earth-centered  axes. 

XAN  =  -sin  (ALAT)  cos  (ALONG)  cos  (ANEL)  cos  (ANAZ)  -sin  (ALONG) 
cos  (ANEL)  sin  (ANAZ)  +  cos  (ALAT)  cos  (ALONG)  sin  (ANEL) 

YAN  =  sin  (ALAT)  sin  (ALONG)  cos  (ANEL)  cos  (ANAZ)  +  cos  (ALONG) 

cos  (ANEL)  sin  (ANAZ)  +  cos  (ALAT)  sin  (ALONG)  sin  (ANEL)  (63) 

ZAN  =  (cos  (ANEL)  (cos  ANEZ)  (cos  ALAT)  +  sin  (ANEL)  sin  ALAT) 
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RESOLUTIONS  FOR  SPIN  STABILIZED  ANTENNA 


The  direction  of  the  spin-stabilized  antenna  is  given  by  the  spin  9'"’=' 
attitude  at  burnout  (BIG)  and  the  spin  axis  heading  at  burnout  (BIH).  B..  tii  of 
these  angles  are  defined  relative  to  the  local  horizontal  plane  with  heading 
defined  positive  clockwise  from  north  and  attitude  defined  positive  in  the  up 
direction.  Two  resolutions  are  required  to  transform  the  unit  spin  axis  vector 
into  components  in  earth-center  coordinates.  One  resolves  into  the  local  tangent 
plane  coordinates  at  the  burnout  condition:  X'  pointing  north;  Y'  pointing  east; 
and  Z'  pointing  up.  These  intermediate  components  are  then  resolved  into 
the  earth-fixed  axis  components,  XI,  YI,  ZI. 


The  axis  system  X'  ,  Y'  ,  Z'  is  fixed  in  inertia  space  and  is  defined  by  the 
burnout  conditions,  i.  e.  ,  the  latitude  at  burnout,  BOLAT,  and  the  longitude  at 
burnout,  BOLONG.  Since  the  earth  revolves,  however,  the  longitude  of  an 
earth-fixed  axis  centered  in  the  earth  relative  to  axis  fixed  in  inertia  space  will 
be  (BOLONG  -  WE  *  TIME)  at  any  given  time  where  (WE)  (TIME),  WET,  is  the 
longitude  change  of  a  point  on  the  earth  due  to  earth  rotation.  The  required 
resolutions  for  the  spin- stabilized  unit  vector  are: 


Resolution  from  unit  vector  e^  to  local  horizontal -plane  (inertial) 
coordinates  X'  ,  Y'  ,  Z'  . 


X' 

cos  BIG  cos  BIH 

"l 

Y' 

= 

cos  BIG  sin  BIH 

0 

Z' 

sin  BIG 

0 

(64) 
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Resolution  from  local  horizontal  axes  components  X'  ,  Y'  ,  Z'  to  earth- 


centered  axis  components  XL,  YI,  ZI. 


XI 

YI 

= 

ZI 

-sin(BOLAT)  cos(BOLONG-WET  -sin(BOLONG-WET)  cos(BOLAT)  cos(BOLONG-WET) 
-sin(BOLAT)  sin(BOLONG-WET)  cos(BOLONG-WET)  cos(BOLAT)  sin(BOLONG-WET) 
cos(BOLAT)  O  sin(BOLAT) 


X 

Y 

Z 


(65) 


XI  =  -sin  (BOLAT)  cos  (BOLONG-WET)  cos  BIG  cos  BIH 
-sin  (BOLONG-WET)  cos  BIG  sin  BIH 
+COS  (BOLAT)  cos  (BOLONG-WET)  sin  BIG 

YI  =  -sin  (BOLAT)  sin  (BOLONG-WET)  cos  BIG  cos  BIH 
+COS  (BOLONG-WET)  cos  BIG  sin  BIH 

+COS  (BOLAT)  sin  (BOLONG-WET)  sin  BIG  (66) 

ZI  ==  cos  (BOLAT)  cos  BIG  cos  BIH  +sin  (BOLAT)  sin  BIG 


If  these  components  are  arranged  to  that  the  time-dependent  terms  are  collected, 
the  expression  becomes 

XI  =  cos  (BOLONG-WET)  [cos  BOLAT  sin  BIG 
-sin  (BOLAT)  cos  BIG  cos  BIH] 

-sin  (BOLONG-WET)  [cos  BIG  sin  BIH] 
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YI  =  cos  (BO LONG-WET)  [cos  BIG  sin  BIH] 

+sin  (BOLONG-WET)  [cos  BOLAT  sin  BIG 

-sin  BOLAT  cos  BIG  cos  BIH]  (67) 

ZI  =  cos  (BOLAT)  cos  BIG  cos  BIH  +sin  (BOLAT)  sin  BIG 

The  terms  enclosed  in  the  brackets  are  constant  depending  only  on  the  burnout 
conditions  BOLAT,  BIG  and  BIH.  They  need  be  computed  only  once  for  any  given 
set  of  initial  conditions. 

RESOLUTIONS  FOR  VEHICLE  ANTENNA  STABILIZED  TO  EARTH  CENTER 

In  this  configuration,  the  vehicle  antenna  is  stabilized  to  point  in  the 
direction  of  the  earth  center  during  the  orbit.  The  resolutions  for  this  case  are 
simply  given  in  terms  of  the  latitude  and  longitude  at  a  particular  time  in  the 
orbit. 


One  resolution  transforms  the  unit  vector  e^^,  into  components  XF,  YF, 
and  ZF  along  the  earth-centered  axis. 


-COS  [LATITUDE  (t.)]  cos  [  LONGITUDE  (t. )] 

1  1 

”®F 

’^F 

= 

-cos  [  LATITUDE  (t. )]  sin  [  LONGITUDE  (t. )] 

1  1 

0 

^F 

-sin  [  LATITUDE  (t.)] 

1 

0 

Where  -e^  is  a  unit  vector  along  a  line  through  the  earth  center  from  the  vehicle 
in  orbit.  Since  the  axis  system  in  the  vehicle  is  defined  with  the  Z-axis  pointing 
upward,  the  e-vector  points  in  the  negative  direction. 
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X_  =  cos  (RLAT)  cos  (RLONG) 

F 

=  cos  (RLAT)  sin  (RLONG)  ^  (69) 

F 

Z_  =  sin  (RLAT) 

F 

RESOLUTION  OF  LINE  OF  SIGHT  FROM  STATION  TO  VEHICLE  (RHO) 

The  direction  of  line  of  sight  to  the  vehicle  is  completely  defined  relative 

to  the  local  tangent  plane  by  the  elevation  angle  (ELE)of  the  line  of  sight 

measured  positive  up  from  the  horizontal  plane  and  by  the  azimuth  or  heading 

angle  (AZR)  measured  positive  clockwise  from  north.  As  before,  two  resolutions 

are  required  to  transform  the  unit  vector  along  the  line  of  sight  into  components 

along  earth-fixed  axes  centered  at  the  earth.  The  first  resolution  transforms 

the  unit  vector  e  into  components  along  the  horizontal  plane  axes  X'  ,  Y'  , 
RHO 

and  Z*  .  The  second  resolution  transforms  this  intermediate  component  into 
components  XRHO,  YRHO,  ZRHO  along  the  earth-fixed  axes  X,  Y,  Z. 

The  first  resolution  is  from  the  unit  vector  e  to  components  X,  Y,  Z 

RHO 

along  local  horizontal -plane  axes. 


X' 

cos  (ELE)  cos  AZR 

®RHO 

Y' 

= 

cos  (ELE)  sin  AZR 

O 

(70) 

Z' 

sin  (ELE) 

O 
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The  second  resolution  is  from  local  horizontal-axis  components  X'  ,  Y'  , 


and  Z'  to  earth-centered  components  XRHO,  YRHO,  and  ZRHO. 


XRHO 

-sin(ALAT)  cos(ALONG)  -sin(ALONG)  +cos(ALAT)  cos(ALONG) 

X' 

YRHO 

= 

-sin(ALAT)  sin(ALONG)  cos(ALONG)  +  cos(ALAT)  sin(ALONG) 

Y' 

ZRHO 

cos  (ALAT)  0  sin(ALAT) 

Z' 

(71) 


or  combining  these  two  matrices  gives 

XRHO  =  -sin  (ALAT)  cos  (ALONG)  cos  (ELE)  cos  AZR 
-sin  (ALONG)  cos  (ELE)  sin  AZR 
+sin  ELE  cos  (ALAT)  cos  (ALONG) 

YRHO  =  -sin  (ALAT)  sin  (ALONG)  cos  (ELE)  cos  AZR 

+COS  (ALONG)  cos  (ELE)  sin  AZR  (72) 

+sin  (ELE)  cos  (ALAT)  sin  (ALONG) 

ZRHO  =  cos  (ALAT)  cos  (ELE)  cos  AZR  +sin  ELE  sin  ALAT 

In  these  components  the  sines  and  cosines  of  ELE  and  AZR  are  the  time-varying 
terms.  ALAT  and  ALONG  are  fixed  by  the  station  location. 
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Finally,  the  cosine  of  the  orientation  angles  between  the  station-  and  vehicle- 
antenna  axes  for  the  various  configurations  can  be  found  by  taking  the  appropriate 
dot  products 


cos  (ALPHA)  =  ■ 

i  =  (XRHO)  (XAN) 

+  (YRHO)  •  (YAN)  +  (ZRHO)  (ZAN) 

(73) 

cos  (ANGLE)  =  ' 

=  (XRHO)  (XI) 

+  (YRHO)  (YI)  +  (ZRHO)  (ZI) 

(74) 

cos  (FNGLE)  =  ‘ 

ip  =(XRHO)(XF) 

+  (YRHO)  (YF)  +  (ZRHO)  (ZF) 

(75) 

Wliere  the  components  of  the  various  unit  vectors  along  the  earth-fixed  axis  are 
given  above  in  terms  of  vehicle  and  station  position. 
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